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Preface 


After getting my Ph.D. in 2014, Icame to Kyoto University for postdoctoral research. 
My host, Prof. Masayoshi Nakashima, invited me to participate in a research project 
on free-standing structures. The free-standing structure has the configuration that 
the superstructure is disconnected from the foundation, and the interface between 
the superstructure base and the foundation is lubricated with graphite to reduce the 
friction coefficient. This idea is very close to sliding isolation. Professor Nakashima 
intended to use this technology in structural retrofitting, so he used a new name. 

After performing a preliminary literature review, I found that scientific research 
on sliding base (SB) structures (I used this name in my later papers, and in some early 
papers, "sliding structures" is commonly used) started as early as the 1980s. Both 
Westermo and Udwadia and Mostaghel et al. published a paper on the response of 
SB structures subjected to harmonic excitation in 1983. Afterward, other researchers 
(e.g., Jangid R. S. at the Indian Institute of Technology) also conducted several analyt- 
ical studies on SB structures. However, there still existed some research gaps in this 
topic prior to 2014. For instance, the theoretical solutions for the sliding-sliding case 
had not been derived; the former studies only used a small number of ground motion 
records, and simplified equations for estimating the peak superstructure response and 
the peak sliding displacement had not been developed. 

During my one year and two months stay in Japan, I completed a study on the 
dynamic responses of two-degree-of-freedom SB systems subjected to harmonic 
excitations. Meanwhile, I took part in the shaking table tests of free-standing 
structures, which was a unique research experience for me. 

After completing my postdoctoral fellowship, I joined Huaqiao University as a 
faculty. My research focus returned to the steel-concrete composite structures, my 
main research interest. However, I still spare some time to continue the left work on 
SB structures. From 2016 to 2022, I, working with two of my graduate students, Fan 
Lin and Li-Wen Xu, completed studies on the peak responses and design methods of 
SB structures subjected to three-component earthquake excitations. My colleague, 
Prof. Yi-Chao Gao, helped me develop the program for response history analyses of 
SB structures subjected to three-component excitations. I am very grateful for his 
help. 


vi Preface 
I wrote this book based on my research on SB structures in the past eight years. 
I hope this book’s publication can guide the design of SB structures and promote 


practical applications of SB structures. 


Xiamen, China Hong-Song Hu 
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Chapter 1 A) 
Introduction EE 


1.1 Fundamentals of Sliding Base (SB) Structures 


Base isolation is an effective approach for reducing damages to structures and their 
contents under severe earthquake excitations. It is generally implemented by using 
special isolators, such as laminated rubber bearings (e.g., Kelly, 1986; Kikuchi & 
Aiken, 1997; Skinner et al., 1993; Yamamoto et al., 2009) and friction pendulum (FP) 
bearings (e.g., Becker & Mahin, 2013; Castaldo & Tubaldi, 2015; Mokha et al., 1991; 
Roussis & Constantinou, 2006). The properties of these isolators can be elaborately 
designed to achieve a certain structural performance. However, they are expensive and 
require high construction techniques. Therefore, when the cost is a major concern, 
base isolation using isolators may not be an appropriate choice. 

Adopting a sliding interface between the base of the superstructure and the foun- 
dation (Fig. 1.1) can also reduce the seismic response of the superstructure. The 
mechanism is very simple: as the friction force between the superstructure and the 
foundation has an upper limit, the seismic force transmitted to the superstructure is 
limited. Structures that adopt this type of isolation technique are called sliding base 
(SB) structures in this book. Since the implementation of SB structures is simple 
and cost effective, they are applicable to low-rise buildings in rural areas. Actually, 
SB structures have been used in some low-rise masonry buildings (Li, 1984; Zhou, 
1997). 


1.2 Practical Implementations of SB Structures 


For the past four decades, several materials have been investigated regarding their 
potential use along the sliding interface of SB structures. Qamaruddin et al. (1986) 
conducted shaking table tests on sliding brick building models with different sliding 
layer materials, namely, graphite powder, dry sand, and wet sand, and obtained fric- 
tion coefficients of 0.25, 0.34 and 0.41, respectively, for the corresponding interfaces. 
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Sliding surface 


Fig. 1.1 Schematic plot of a building adopting SB technique 


Tehrani and Hasani (1996) conducted experimental studies on adobe buildings with 
dune sand and lightweight expanded clay as sliding layers; the friction coefficients 
were reported as 0.25 for dune sand and 0.2—0.3 for lightweight expanded clay. Barba- 
gallo et al. (2017) tested a steel-mortar interface lubricated by graphite powder; the 
static and dynamic friction coefficients were close to 0.19 and 0.16, respectively, and 
they were independent of both the sliding velocity and the superstructure properties. 

Polymer materials are also suitable choices for the sliding interface. Yegian 
et al. (2004) investigated the frictional characteristics of four synthetic interfaces 
[namely, geotextile-high density polyethylene (HDPE), polypropylene (PTFE)- 
PTFE, ultrahigh molecular weight polyethylene (UHMWPE)-UHMWPE, and 
geotextile-UHMWPE] as potential candidates for sliding isolation through cyclic 
and shaking table tests. It was determined that the geotextile-UHMWPE interface 
was suitable for sliding isolation applications because the friction coefficient of this 
interface is insensitive to large variations in the sliding velocity and normal stress; as 
a result, this interface can easily be introduced into engineering design. The obtained 
static and dynamic friction coefficients of the geotextile-UHMWPE interface were 
approximately 0.11 and 0.08, respectively. Nanda et al. (2012, 2015) conducted 
experimental studies on four sliding interfaces with green marble against HDPE, 
green marble, geosynthetics and rubber sheeting, respectively. The static friction 
coefficients were found to be independent of the normal stress, and the dynamic 
friction coefficients were insensitive to variations in the sliding velocity. Moreover, 
the observed dynamic friction coefficients of the four investigated interfaces ranged 
from 0.07 to 0.18, and the relative differences between the static and dynamic friction 
coefficients were all below 15%. Jampole et al. (2016) adopted sliding isolation bear- 
ings consisting of HDPE sliders and galvanized steel surfaces to seismically isolate 
light-frame residential houses; shaking table tests showed that the friction coefficient 
of this sliding interface was nearly 0.18 with a slight variation between the stick and 
sliding phases. 
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The aforementioned sliding interfaces were basically insensitive to the variations 
in the sliding velocity and pressure. Therefore, the Coulomb friction model can be 
used to model the behavior of the sliding interface of an SB structure. 


13 Review of Analytical Studies on SB Structures 


1.3.1 Studies on 2DOF SB Systems Subjected to Harmonic 
Ground Motions 


The simplest model for an SB structure contains two masses, one for the superstruc- 
ture and the other for the sliding base; thus, this model is referred to as a 2-degree- 
of-freedom (2DOF) SB system. Westermo and Udwadia (1983) and Mostaghel et al. 
(1983) first studied the dynamic responses of 2DOF SB systems under harmonic exci- 
tations, and both groups developed the governing equations of motion and numer- 
ical implementations needed to perform response history analyses of such systems. 
Westermo and Udwadia (1983) pointed out that the response of this system under 
harmonic excitations converged rapidly to a periodic response after several cycles. 
Three different periodic responses were observed, namely stick-stick, stick-sliding, 
and sliding-sliding cases, depending on the amplitude of the input accelerations 
and the structural characteristics of the system. They also derived the explicit equa- 
tions for the condition of the initiation of the stick-sliding case. Mostaghel et al. 
(1983) conducted parametric studies for the critical responses of the 2DOF SB 
system under harmonic ground motions. Iura et al. (1992) followed the work of West- 
ermo and Udwadia (1983), and Mostaghel et al. (1983). They derived the analytical 
expressions for the condition of the initiation of the sliding-sliding case. Therefore, 
combined with the work of Westermo and Udwadia (1983), the explicit equations 
for the occurrence conditions of three periodic response cases were obtained. More 
recently, Hu and Nakashima (2017) conducted a comprehensive parametric study on 
the maximum responses of 2DOF SB systems under harmonic ground motions and 
derived a theoretical solution for the response corresponding to the sliding-sliding 
case. 


1.3.2 Studies on SB Structures Under Earthquake Excitation 


Mostaghel and Tanbakuchi (1983) studied the seismic responses of 2DOF SB systems 
subjected to the N-S component of the El-Centro record from the 1940 Imperial 
Valley earthquake and the S86E component of the Olympia record from the 1949 
Western Washington earthquake. Response spectra of the absolute acceleration and 
sliding displacement were developed. It was found that the SB isolation can effec- 
tively control the level of the superstructure response. Qamaruddin et al. (1986b) 
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conducted a study similar to that of Mostaghel and Tanbakuchi (1983), with emphasis 
on using the SB system in masonry buildings. The N-S component of the El Centro 
record and the longitudinal component of the Koyna record were considered. The 
findings of their study were similar to those of Mostaghel and Tanbakuchi (1983). 
Yang et al. (1990) and Vafai et al. (2001) studied the responses of multiple-degree- 
of-freedom (MDOF) SB structures. Their focus was to develop efficient numerical 
methods for response history analyses, whereas the response characteristics of these 
structures were not sufficiently addressed. 

Jangid (1996a) compared the responses of single-story SB structures subjected 
to two components and a single component of the El-Centro record; the numer- 
ical results indicated that the former increased the sliding displacement and reduced 
the absolute acceleration of the superstructure in comparison with the latter. Shakib 
and Fuladgar (2003a) adopted the same model as Jangid (19962) but included the 
vertical component in the ground motion input, and three ground motion records 
were considered, namely, the El-Centro record, the Tabas record from the 1978 
Tabas earthquake, and the Renaldi record from the 1994 Northridge earthquake. The 
effect of the vertical component was highly dependent on the superstructure period 
and the input ground motions; additionally, the responses of low-period structures 
could be strongly affected by the vertical component of the ground motion, while 
this influence was insignificant when the superstructure period exceeded 0.7 s. Jangid 
(1996b) and Shakib and Fuladgar (2003b) also studied the responses of asymmetric 
single-story SB structures under multidimensional inputs. Both studies indicated that 
the bidirectional interaction between the frictional resistance at the sliding interface 
and the vertical component of the earthquake excitation could significantly affect the 
responses of torsionally coupled systems with SB isolation. More recently, Hu et al. 
(2020, 2022) investigated the peak superstructure response and peak sliding displace- 
ment of SB structures subjected to three-component earthquake excitation using 
a large number of ground motion records. The influence of various structural and 
ground motion characteristics on these two response quantities was comprehensively 
studied, and simplified design equations were also developed. 
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Chapter 2 A) 
Responses of 2DOF Sliding Base Systems |as 
Under Harmonic Ground Motions 


2.4 Equations of Motion 


For a 2DOF SB system shown in Fig. 2.1, the dynamic equilibrium equations can be 
written as 


DE (2.1) 
Mp(iig + tis) + m(iig + üs + ii) = f 
in which m, k, and c refer to the top mass, lateral stiffness, and viscous damping 
coefficient of the superstructure, respectively; mp is the mass of the sliding base; 
Ug(t), Us(t), and u,(t) are the ground displacement, sliding displacement, and relative 
displacement between the top mass and sliding base, respectively; i, us, and ú, 
are the corresponding velocities; 4, is, and ü, are the corresponding; and f is the 
friction force between the sliding base and foundation. The first equation of Eq. (2.1) 
denotes the dynamic equilibrium of the top mass, while the second equation is the 
dynamic equilibrium of the entire system. 

The SB system can display two kinds of phases in its response history: the stick 
phase and the sliding phase. For the stick phases, the sliding acceleration, iis, is 
equal to 0, and the sliding friction force is greater than the friction force, f; therefore, 
Eq. (2.1) leads to 


mu, + cu, + ku, = —mü, 
aii, + üg = f/(m + my) (2.2) 
|f] < (m + mp) ug 
where 
ui (2.3) 
a= : 
m+mp 
© The Author(s) 2023 5 
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Fig. 2.1 A 2DOF sliding Top mass, m 
base system A 


is the mass ratio; yz is the friction coefficient; and g is the gravity acceleration. The 
first equation of Eq. (2.2) can be written as 


ü, + 26wu, + wu, = —ii, (2.4) 
where 
k 
pa e (2.5) 
m 2am 


o and & are the natural frequency and damping ratio of the corresponding fixed 
base structure, respectively. Equation (2.4) is the differential equation governing the 
relative displacement of a single-degree-of-freedom (SDOF) system under ground 
acceleration, ig. 

According to the second and third equations of Eq. (2.2), we obtain 


aü, -- ü,| < ug (2.6) 


This is the precondition for the stick phases. Sliding occurs when it is no longer 
satisfied, and 


f —ó(m-my)ug (2.7) 


in which 6 represents the direction of the friction force, and |ó| = 1. While transi- 
tioning from the stick phase to the sliding phase, the direction of the friction force 
remains unchanged. Thus, it can be inferred from the second equation of Eq. (2.2) 
that ô and wii, + ij, have the same sign at the moment before sliding. By substituting 


2.1 Equations of Motion 7 
Eq. (2.7) into the second equation of Eq. (2.1), we obtain 
li, = dug — ad, — üg (2.8) 
Substituting Eq. (2.8) into the first equation of Eq. (2.1) to obtain 
(1 —a)mu, + cu, + ku, = —óumg (2.9) 


Dividing Eq. (2.9) by (1 — a)m to obtain 


ô 
ü, + 26e, + ofu = -5 (2.10) 
—-a 
where 
w £ 
ei = = (2.11) 
l-a l-a 


Equation (2.10) is the differential equation of a SDOF system with the natural 
frequency of w; and damping ratio of 4; under a step force. The static displacement 
corresponding to the step force is 


ig = E (2.12) 


Thus, in the sliding phases, a new natural frequency, o, and a new damping ratio, 
£i, of the relative displacement vibration are determined, both of which are linked to 
the mass ratio, a. The solution of Eq. (2.10) is 


u, (t) = u, (ti) 
Coy T aof ue) TO) (2.13) 


in which ft is the global time, t; is the moment when sliding starts, t = t — t;, and 
A(T) = e hit 


2 sin(wiaT) 

cos(@jgt) + E1/4/1 — ££ sin(@jaT EL 
(iat) + &1/y/1 — Et sinoat) "IE 
iom cos(wrat) — &1/y/ 1 — Ep sin(onar) 

791 


(2.14) 


1— eer C 1 — & sin(@yat) + eso) 
b(t) = us; E 
e 5e o sin(wiaT) 
1—4? 


(2.15) 
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Yes No 
ga No Determing the m 
Non-sliding ii p> sliding direction —> ą Sliding 


Yes 


Fig. 2.2 Flow chart for calculating the response history of the 2DOF sliding base system 


wid = wiy 1 — &? (2.16) 


The first term in the right-hand side of Eq. (2.13) corresponds to free vibration 
caused by the initial condition, while the second term represents the vibration caused 
by the step force. 

At the onset of sliding, the velocity of sliding equals to 0, i.e., us(t;) = 0. 
Therefore, integrating Eq. (2.8) leads to 


ùs (t) = dug(t — ti) — a (Ùr (t) — ùr (ti) — (g(t) — g(t) (2.17) 


this round Sliding stops once 4, (1) returns to 0 again. The sliding of the structure can 
either stop or persist, based on whether Eq. (2.6) is satisfied or not. For this check of 
Eq. (2.6), Eq. (2.4) should be used to determine the relative acceleration ij,, under 
the assumption that sliding stops. The sliding displacement, u; (ft), can be obtained 
by integrating Eq. (2.17): 


Us(t) = us(ti) + A — t)! — o(u, (t) — ur (ti) — à, (t) — t;) 
(ug(t) — ug (ti) — ùg (ti) — ti)) (2.18) 


Figure 2.2 summarizes the process of computing the response history. 


2.2 Typical Response Histories 


In the following analyses, the ground acceleration is taken as a sinusoidal wave if 
not specified: 


ii, = a, Sin(wet) (2.19) 


where a, and c, are the amplitude and frequency of the sinusoidal wave, respectively. 
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Figure 2.3 shows the responses of a sliding base system with u = 0.2, a = 0.5, 
T — 2x [o =0.5 s, E =5%,T, = 2T /w, = ls, and a, = 0.4g, in which T represents 
the natural period of the superstructure, and T, indicates the period of the ground 
acceleration. As shown in Fig. 2.3, the responses of the system converge to steady 
periodic responses with the same period as the ground acceleration after several 
cycles. The sliding base system pauses briefly before sliding in the opposite direction 
in this instance. This is referred to as the stick-sliding case. After the amplitude of the 
ground acceleration exceeds a certain value, the sliding base system will continue 
to slide incessantly during the steady periodic state, as shown in Fig. 2.4, where the 
ground acceleration amplitude, ag, is increased to 1.2g. This is referred to as the 
sliding-sliding case. Sliding will not occur if the ground acceleration is sufficiently 
low, which is known as stick-stick case. The sliding base structure and the fixed base 
structure have no difference in this case. The following analyses focus on the steady 
periodic responses. 
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Fig. 2.3 Responses of the sliding base system for u = 0.2, a = 0.5, T = 2m /w = 0.5 s, & = 5%, 
T, =21/wg = 1 s, a, = 0.4g: a absolute acceleration of the top mass; and b sliding displacement 
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Fig. 2.4 Responses of the sliding base system for u = 0.2, a 
T, —2z/og = 1 S, dg 


0.5, T = 2x /w = 0.5 s, £ = 5%, 
1.2g: a absolute acceleration of the top mass; and b sliding displacement 


2.3 Occurrence Conditions of the Three Types of Periodic 
Responses 


2.3.1 Boundaries Between the Stick-Stick and Stick-Sliding 
Cases 


For the stick-stick case, the steady response of the relative displacement is 


75 Ru sin(wgt — $) (2.20) 


ae (2:21) 
t= (o/o) 


where ¢ is the phase angle, which defines the time by which the response lags behind 


the input ground motion. Substituting Eqs. (2.19) and (2.20) into the precondition 
for the stick phases (i.e., Eq. 2.6) gives 


lo R, sin(wet — p) + sin(wet)| < wg/ag (2.22) 


2.3 Occurrence Conditions of the Three Types of Periodic Responses 11 
where 


2 
Ry = (w/w) Ra = (e/o) (2.23) 


JIi- oro + Beloof 


is the acceleration response factor for the fixed base structures. Equation (2.22) can 
be rewritten as 


|(o Ra cos $ + 1) sin(wet) — æ Ra sin $ cos(wet)| < ug/as (2.24) 


Equation (2.24) is always satisfied for the stick-stick cases, so the maximum value 
of the left-hand side term should be smaller than jxg/a,; thus, we have 


Ag/dg > V (eR, +2(aR,) cos + 1 (2.25) 
Therefore, the boundaries between the stick-stick and stick-sliding cases are 


dg 1 


ug V(a@Ra)? + 2(aR,) coso + 1 


(2.26) 


The variations of the boundaries between the stick-stick and stick-sliding cases 
for different w/w and o are shown in Fig. 2.5. 


Fig. 2.5 Boundaries 30 
between the stick-stick and 
stick-sliding cases (£ = 5%) 25 
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2.3.2 Boundaries Between the Stick-Sliding 
and Sliding-Sliding Cases 


As derived in Sect. 2.2, Eq. (2.10) is the governing equation of the relative displace- 
ment during the sliding phase. Equation (2.10) is the differential equation of a SDOF 
fixed base system with the natural frequency of w; and damping ratio of 5, subjected 
to a step force. The static displacement corresponding to this step force, tsr, is given 
in Eq. (2.12), which has the same sign as the sliding direction. As shown in Fig. 2.4, 
for the sliding-sliding case, the sliding base system slides in one direction for half 
period of the ground motion, 0.57,, and then slides in the opposite direction for 
another 0.5T',, and so on. Therefore, the equivalent step force also changes direction 
in 0.57, when the sliding direction changes, as shown in Fig. 2.6. Thus, the relative 
displacements of the two opposite half sliding cycles have the same absolute values 
but opposite signs. 

For a certain half sliding cycle, the solution of Eq. (2.10) (i.e. Eq. 2.13) can be 


rewritten as 
e = ao(s) + b(t) (2.27) 
u, (t) uro 


where u,o and 1,9 are the relative displacement and velocity at the moment when this 
half sliding cycle starts, and this moment is denoted as f;; T = t — t; is the local time 
during this half sliding cycle, so O € t < m/w,. When t = 7 /@g, the opposite half 
sliding cycle starts, the relative displacement and velocity at this moment are equal 
to —u,o and —u,o, respectively. Therefore, 


u(t) 
Plk: PO mC 
— 


057, 05% 057 


I KI» I Jemma 


-ugla 


Fig. 2.6 Dynamic model for the relative displacement 


2.3 Occurrence Conditions of the Three Types of Periodic Responses 


(S2) = aaroo ( 12) +o) 


The solutions of Eq. (2.28) are 


— sinh 0; + (s Lib s?) sin 6 


cosh 6; + cos 62 


nu sinh 0, 4- (sv t= s) sin 62 


o»? cosh 0; + cos 05 


Uro = Ust 


sin 05/ 1-& Suga sin 05 / 1-& 


lo = —UstW, = 2 
cosh 6; + cos 62 œ coshdé; + cos 65 
where 
[2 [2 Mid 
0; = §— 7, h =,/1- 2— r = — r 
Wg [5 Og 


(2.28) 


(2.29) 


(2.30) 


(2.31) 


and ô represents the direction of the friction force, i.e., the opposite direction of 


sliding. 


By substituting Eqs. (2.29) and (2.30) into Eq. (2.27), the relative displacement 


and velocity at any time of this half sliding cycle are obtained: 
ur +2) = E 
@ 


1—e 81 @1t 


cos(05 — wjgTt),/ 1 E? &, sin(05 — wjgt) + & sin(@jqt)e®! +,/1-— E? cos(@jqt)e®! 


1— E? (cosh 01 + cos 05) 


-ô sin(c41) (e^ + cos 02) — cos(w1aT) sin 0 
ü, (ti +T) = HOY Eoi (1a X ») (waT) 2 


2 
o 1 — £? (cosh 0; + cos 62) 


have 


üs (t; + I/We) = óngr/o, — a(— 21,9) — 2ag/W, cos(wgti) =0 


(2.32) 


(2.33) 


When t = 7 /@g, the opposite half sliding cycle initiates. From Eq. (2.17), we 


(2.34) 
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Therefore, 

OT US | AWeltro óng 

cos(w,t;) = + — = Áı (2.35) 
(efi) 
2a, ag ag 
where 

sin 0 T 

A, = LA EEEE (2.36) 


o1- E? cosh; +cosh 2 
y 1 


Since a certain half sliding cycle initiates at t = t;, 
ali, + 77 NN = 8(a(— 2£oi,9 — Xu, — ag sin(@gti)) + dg sin(w,t;)) 
= —A»sug + ó(1 — o)a; sin(oyt;) > ug (2.37) 
where 


a(s sin 62/,/1 — £2 + sinh a) 


Ao = 2.38 
cosh 6; + cos 65 ( ) 


From Eqs. (2.35) and (2.37), we have 
(1 — o)" (a — Aj(ug)?) = + A2 (ug (2.39) 


Thus, the boundaries between the stick-sliding and sliding-sliding cases are 


Ae cU A? / a? + AF (2.40) 


The variations of the boundaries between the stick-sliding and sliding-sliding 
cases for different c, /c and a are shown in Fig. 2.7. 

When c/c, — 0,0; and 0; > +00, A, — 7/2, and A2 — a; thus, Eq. (2.40) 
becomes 


a = V 1 -0!/0 — ay + 0/2? (2.41) 


When w/w — +00, 6; and 0) — 0, A; > zx/(2(1 — a)), and A2 — 0; thus, 
Eq. (2.40) becomes 


ag; vlc (x/2y 


T ERE ur (2.42) 


2.4 Parametric Study for the Maximum Responses 15 


Fig. 2.7 Boundaries 50 
between the stick-sliding and 
sliding-sliding cases (£ = 
5%) 40 
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As can be seen in Fig. 2.7, there are also some local peak values when 62 = (2n — 
1)z. These points are corresponding to the resonance. 


2.4 Parametric Study for the Maximum Responses 


2.4.1 Maximum Pseudo Acceleration of the Top Mass 


The maximum shear force applied to the superstructure can be represented by the 
maximum pseudo acceleration of the top mass (Chopra, 2001), A, as follows: 


A = o? x (max|u,(t)|) = E (maxlu, (D) (2.43) 


The function of dg, 4g, @,/w, a, and £ is related to A in the 2DOF sliding base 
system depicted in Fig. 2.1. Figure 2.8 illustrates how the normalized maximum 
pseudo acceleration, A/ug, relates to the frequency ratio, c/c, for varying values 
of a and ag/jg. A sliding base structure, unlike fixed base structures, has multiple 
resonant frequencies when sliding occurs. a, / ug and o can also affect these resonant 
frequencies. As shown in Fig. 2.8, b, as a, / ug increases, several resonant frequency 
ratios appear in the region of w,/w < 1, and these ratios shift towards higher values 
gradually and eventually reach upper limits, which correspond to the sliding-sliding 
cases. As shown in Fig. 2.8c, d, resonances are more prone to happen with smaller 
mass ratios, o. 

As shown in Fig. 2.8c, d, the maximum pseudo acceleration typically reduces 
with an increase in the mass ratio, o. The reason for this is that the damping ratio 
corresponding to the sliding phase, &), is equal to £/4/1 — o as given in Eq. (2.11), 
which increases as o increases. However, in the resonant frequency ranges of a 
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Fig. 2.8 Relationship between the normalized maximum pseudo acceleration and wg /c (& = 5%): 
a œ —05;bo = 0.8; € ag / ug = 2; and d a,/ug = 8 


specific sliding base structure, the maximum pseudo acceleration may be greater 
compared to other structures with lower mass ratios. 

Figure 2.9 shows the relationship between A/ug and a,/jg for different values 
of a and c/c, where the circle and triangle represent the boundary between the 
stick-stick and stick-sliding cases, and the boundary between the stick-sliding and 
sliding-sliding cases, respectively. The A/jg versus a,/j1g curves are segregated 
into three areas, each corresponding to cases of stick-stick, stick-sliding, and sliding- 
sliding. In the stick-stick region, responses remain independent of mass ratio, o, with 
the A/ug versus a/g curve taking the form of a straight line. The slope of this 
straight line is equal to the displacement response factor, Ry, for a fixed base structure: 


1 
Ry = (2.44) 


IE es] + pitur 


In the stick-sliding region, due to additional resonances when w,/w < 1, the 
maximum pseudo acceleration of a sliding base structure may exceed that of the 
related fixed base structure. This phenomenon can be observed in Fig. 2.9a for w,/w 
= 0.1 and 0.5. In the sliding-sliding region, the upper limit of the maximum pseudo 


2.4 Parametric Study for the Maximum Responses 17 


(a) asl ug 


O Boundary between the stick-stick and stick-sliding cases 


^ Boundary between the stick-sliding and sliding-sliding cases 


Fig. 2.9 Relationship between the normalized maximum pseudo acceleration and a; / ug (E = 596): 
aa = 0.5; b w/w = 0.5 


acceleration is reached, and does not change even when there are larger ground 
accelerations. This upper limit reflects the effectiveness of sliding base structures 
for isolating extremely large earthquakes. As shown in Fig. 2.9b, in normal circum- 
stances, as the mass ratio increases, the likelihood of the sliding-sliding case occurring 
becomes more difficult. This trend can be more clearly observed in Fig. 2.7. 


2.4. Amplitude of the Sliding Displacement 


As shown in Fig. 2.4b, the maximum sliding displacement results from the accu- 
mulation of transient responses that occur before reaching the steady periodic state. 
Hence, the maximum sliding displacement is considerably influenced by the initial 
ground acceleration, e.g., there is a significant disparity in the maximum sliding 
displacement between the sine and cosine ground accelerations. To represent the 
extent of sliding when exposed to harmonic ground motions, the amplitude of the 
sliding displacement, u;,), is a more appropriate response quantity compared to the 
maximum sliding displacement. It is defined as the difference between the maximum 
and minimum sliding displacements during the steady state; therefore, the value of 
it is exclusively linked to the responses of the steady state. 

Figure 2.10 depicts the correlation between the frequency ratio, c/c, and the 
normalized sliding displacement amplitude, ws ap / (ag / 0»), for varying a and a,/ug 
values. As shown in Fig. 2.10a, b, when o /c A is less than 1, the normalized sliding 
displacement amplitude is not greatly affected by the mass, œ, ratio or frequency ratio, 
€ / c. There is a noticeable reduction of the sliding displacement amplitude when 
€ / increases, after the frequency ratio, c /c, surpasses 1. This result is in agree- 
ment with Fig. 2.5, where the critical a,/1g for the occurrence of the stick-sliding 
case increases as w, /w increases once c; /w surpasses 1. As shown in Fig. 2.10b, if 
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Fig. 2.10 Relationship between the normalized sliding displacement amplitude and wg/w (E = 
590):aag/ug = 2; and b ag /ug = 8 


the normalized ground acceleration amplitude, a, / ug, is sufficiently large, as w,/w 
surpasses a certain value, the sliding displacement amplitude will rise alongside an 
increase in c/o. Figure 2.5 is also in accordance with this, where the critical value 
of ag/jg leading to stick-sliding decreases with increasing value of c, /c once the 
peak point of critical value a, / ug is reached. Figure 2.10b also indicates that there is 
a tendency generally for the sliding displacement amplitude to decrease as the mass 
ratio increases in the area of large frequency ratios. 

Figure 2.11 depicts the correlation between the normalized sliding displacement 
amplitude, us ap/ (a mi w?) , and the normalized ground acceleration amplitude, a, / ug, 
for different w; /c and a. As shown in Fig. 2.11, after sliding occurs, the normalized 
sliding displacement amplitude increases as a,/jg increases, but with a gradually 
decreasing speed increase. The normalized sliding displacement amplitude has an 
upper bound value that is 2. The reason for this is that when a,/g is very large, 
the absolute acceleration of the sliding base (which is limited by the friction coeffi- 
cient) is negligible compared with the ground acceleration, so the sliding base can 
be considered motionless from the perspective of ground, and us ap is equal to the 
vibration amplitude of the ground displacement, 2a, / o». 


2.5 "Theoretical Solutions for the Responses 
of the Sliding-Sliding Case 


For the sliding-sliding case, the maximum pseudo acceleration reaches the upper 
limit, and is no longer dependent on the amplitude of the ground acceleration. This 
upper limit response has important physical meaning, and reflects the effectiveness 
of the sliding base system for reducing the superstructure response. 

In this section, the theoretical solutions for the maximum pseudo acceleration 
and sliding displacement amplitude for the sliding-sliding case are derived. The 
derived results are further used to explain the mechanism of the sliding base system. 
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Fig. 2.11 Relationship between the normalized sliding displacement amplitude and a/g (§ = 
5%): aa = 0.5; b wg /w — 0.5 


For the stick-sliding case, theoretical solutions cannot be obtained because of the 
complicated transformation between stick and sliding phases. 


2.5.1 Solutions for the Maximum Pseudo Acceleration 


The maximum relative displacement will be obtained when the relative velocity is 
equal to 0, i.e. ú, = 0; thus, the corresponding local time, tj, can be obtained by 
taking the right-hand side of Eq. (2.33) equal to zero: 


sin 05 


tan(wiatj) = e^: + cos 0; 


(2.45) 
From Eq. (2.45), several solutions for t; may be obtained, and the number of 


solutions is dependent on the value of 62. By substituting Eq. (2.45) into Eq. (2.32), 
the relative displacement at the moment of v; can be obtained: 


E 
—óng T uv cos(6» — @\dT;) + cos(wiatj)e"" 


w? cosh 6; + cos 65 


u,(t; + tj) = 


(2.46) 
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Using Eq. (2.46), the normalized maximum pseudo acceleration can be written 
as 


U, (ti + Tj)@ 
HE 


A/ug = max 


E 
is vor cos(@2 — v47;) + cos(wiatj)e" 


(2.47) 
cosh 6; + cos 6; 


= max 


Equation (2.47) is not an explicit equation yet, since there may be more than one 
solution for t; and the corresponding pseudo acceleration. The maximum pseudo 
acceleration is the maximum of these candidates. To obtain the explicit form of 
Eq. (2.47), different cases of 0» need to be considered. The following presents the 
complete derivations for the explicit expressions of the normalized maximum pseudo 
acceleration, A / ug. 

Equation (2.45) leads to 


tan? (wiat) sin? 6) 


= 2.48 
tan?(ci47;j) +1 2e% (cosh 0 + cos 62) oe) 


sin? (wiaty) = 


If 0 A nr (n=1,2,...), sin@2 z 0, so sin(wiatj) Æ 0 from Eq. (2.48). 
Therefore, Eq. (2.47) can be rewritten as Eq. (2.49) by using Eq. (2.45). 


gi 


—— ingo 
A/ug = max l-e Vi | — 2 (2.49) 
sin (via tj)(cosh 01 + cos 2) 
Case 1. If 0 < 62 < x, sin 62 > 0. From Eq. (2.45), we have 
in 0. in 0 
tanfor) = ———— m2 — tan(02/2) (2.50) 


e? --cos0; 1+ cos05 
Since 0 < v; < z/o,,0 < wiaTj X 65, so there is only one solution for tj, and 


sin 05 


—— ——— 05/2 2 2.51 
avez). ee on 


@M1dTj = arctan( 


From Eq. (2.51), we have sin(@1aT;) > 0, so from Eq. (2.48), we have 


in 
Sn? — 269 (cosh 6, + cos6;) (2.52) 
sin(@1aT;) 
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Substituting Eq. (2.52) into Eq. (2.49) gives 


=i =01dTj 2e 
A =|l-e viii — 2.53 
/u8 Í cosh 6; + cos 65 ( ) 


Case 2. If zt < 05 < 27, sin 0 < 0, From Eq. (2.45), we have 


tan(62/2) sin 05 i ( ) sin 05 (2.54) 
an = — — «tan(ej4T;) = ———_ : 
2 1 + cos 0; P e. + cos 05 
so the smallest solution for v; is 
J8 < 2/2 t sé Y. (2.55) 
mw/2< < @ qT; = arctan| ————— m mc : 
" mg e?! + cos 05 
Since 
sin 62 
arctan| —————— ] + 27 > 65/24 1 > 6; (2.56) 
e? + cos 05 


there is no other solutions for v; because 0 < «47; < 62. From Eq. (2.55), we have 
sin(@14T;) > 0, so from Eq. (2.48), we have 


sin 05 


sin(@14T;) 


= —,/2¢% (cosh 6; + cos 62) (2.57) 


Substituting Eq. (2.57) into Eq. (2.49) gives 


E 


-—H. ew; 2e 
A E NIS A 2.58 
/u8 =n cosh 0; + cos 05 ( ) 


Case 3. If 0; > 27 and 0; Æ nz (n = 1,2,...), there are at least two solutions 
for t;. If sin; > 0, i.e., 2nz < 0; < Qn + 1)z, then from Eq. (2.45), we have 
tan(@gT i) > 0. Therefore, the first solution for v; is 


0 t ay /2 (2.59) 
< T; = arctan| ————— <7 » 
de. e?! + cos 05 


and the corresponding pseudo acceleration is 


—~—L¢ 2e 
"C ---———— (2.60) 
cosh 6; + cos 05 


u, (ti + Tj)@” 


HE 
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where 


$ -— arctan ( =") (2.61) 


e?! + cos 0; 
The second solution for r; is 


sin 05 


JT < Ojqtj = arctan ( 5755 


) ta < 3/2 (2.62) 


and the corresponding pseudo acceleration is 


- - 5 +n) efi 
ES! AD — 2.63 
Te cosh 6; + cos 6; ( ) 


The pseudo accelerations corresponding to other solutions of t; are all smaller 
than the larger one of Eqs. (2.60) and (2.63), since the oscillation of the relative 
displacement gradually decays because of the damping. If Eq. (2.60) is equal to 


-—j¢ 2e 
"UT Va LL (2.64) 
cosh 6; + cos 05 


it is clear that Eq. (2.63) is larger than Eq. (2.60). If Eq. (2.60) is equal to 


uy (ti + tj)o? 
HE 


U, (ti + tj)o? 
HE 


- Ke 2e% 
E e 1 (2.65) 
cosh 6; + cos 6; 


Ur (ti + tj)o? 
HE 


Equation (2.60) minus Eq. (2.63) is 
EE! PANEG- ge 
" ae 2^ l-e v \-2 
cosh 6; + cos 6» 


2e% Lg 
< (1 eR") 2 
cosh 6; + cos 6; 
2e% 
| 1—e9)—-2 
= cosh 6; — 1 ( i ) 


1—e*94)-220 (2.66) 


Therefore, Eq. (2.63) is always larger than Eq. (2.60), and the maximum pseudo 
acceleration is obtained at 
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sin 05 


T < @ qT; = arctan| ——— —— 
7 e^ + cos 05 


) ta < 3/2 (2.67) 


and the maximum pseudo acceleration is 


ce =01dTj 2e 
A =1 NIU | —— 2.68 
d TE cosh 0, + cos 63 ( ) 


Case 4. If sin@. < 0, i.e., (2n + 1)a < 0; < (2n + 2)z, then from Eq. (2.45), we 
have tan (wia T i) < 0. Therefore, the maximum pseudo acceleration will be obtained 
at the first solution of t}, that is 


in 6 
7/2 < e41j = arctan mE ta +a <7 (2.69) 
e? + cos 05 
and the maximum pseudo acceleration is 
-— euni 2e: 
A/ug=1te Vv" ——————— (2.70) 


cosh 6; + cos 65 


Case 5. If 0 = nz (n=1,2,...), sinb? = O0, and from Eq. (2.45), we have 
tan(ciat;) = 0, so the first solution for tT; is 


ciat; =0 (2.71) 


Substituting Eq. (2.71) into Eq. (2.46) gives the corresponding pseudo accelera- 
tion: 


u, (t; + Tti) —1)"+e% 
it yor! Dte oi (2.72) 
ug cosh 6, + (— 1) 
The second solution for r; is 
WidTj ST (2.73) 


and Substituting Eq. (2.73) into Eq. (2.46) gives the corresponding pseudo acceler- 
ation: 


ES 


Hite Vo Cute 
B cosh 6; + (— 1)" 


meee ean VY" +e (2.74) 
cosh 6; + (— 1)” 


Ur (ti + Tj)@ 
HE 
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Equation (2.72) minus Eq. (2.74) is 


Kee Oy aim De ae tg 
cosh 6; + (— 1)" ~ cosh 6; + (— 1)" 
—-2 ee : 
_ aie) < 0 (niseven) (2.75) 
0 (n is odd) 


Therefore, Eq. (2.72) is always smaller than Eq. (2.74) For 00 = 
nz (n = 1,2,...), 


(— 1)" + e& 2e 2e% 
n — ni (2.76) 
cosh 0; 4- (— 1) cosh 0; + (— 1) cosh 0; + cos 6; 
so the maximum pseudo acceleration can be written as 
-— r 2e 
A/ug=1+e v (2.77) 


cosh 0; + cos 05 


The results of Cases 2—5 can be combined in one equation, that is, for 62 > 7r, 
the normalized maximum pseudo acceleration is 


-— ouat 2e% 
A ES EH ——————— 2.78 
/H8 cosh 6; + cos 65 ( ) 
where 
in 0 
m/2 < @dTj = arctan ( =) ta <3n/2 (2.79) 


In summary, the explicit solutions for the maximum pseudo acceleration are as 
follows: 
For 0 < 6) < x, the normalized maximum pseudo acceleration is 


= - z O14 To 2e 
A —1- nil ——————— 
d à cosh 6; + cos 65» 


sin 05 ) T, 


(2.80) 


where 


(2.81) 


QjqTt9 = arctan| —————— 
Heg (s + cos 05 
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For 62 > x, the normalized maximum pseudo acceleration is 


A TOV ME a 2.82 
= | + il . 
/H8 E cosh 6; + cos 62 ( ) 


sin 05 
e? + cos 05 


where 


7/2 < 0410 = arctan( ) +a < 3/2 (2.83) 


In the above equations, ro is the local time corresponding to the maximum pseudo 
acceleration. 0, and 62 in Eqs. (2.80)-(2.83) are only related to & and w/w, as 
given in Eq. (2.31); thus, the normalized maximum pseudo acceleration, A /ug, is 
only dependent on á; and wı /c, for the sliding-sliding case. € and @ are the natural 
frequency and damping ratio for the sliding phase, which are related to the mass ratio, 
a, as shown in Eq. (2.11). 


2.5.2 Interpretation of the Solutions for the Maximum 
Pseudo Acceleration 


Figure 2.12 plots the relationship between A/j1g and w,/w using the derived theo- 
retical solutions. As shown in Fig. 2.12a, the results from the theoretical solutions 
and numerical methods are the same, which verifies the accuracy of the theoretical 
solutions. As shown in Fig. 2.12b, as the mass ratio, a, increases, the frequency ratios, 
€ / €, of resonance shift towards larger values and the general responses decrease. 
This is because é; and œ; increase as o increases for given & and w. 


— Theoretical solutions 


=0. 
=0. 
O Numerical results zy 
=0. 


0.1 1 10 
(a) Og CD 


Fig. 2.12 Relationship between A/ug and w/w (§ = 5%): a verification of the theoretical 
solutions for æ = 0.5; b results for different values of a 
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In the following paragraphs, we will examine various frequency ratio regions, 
and elucidate the outcomes illustrated in Fig. 2.12b. When w,/w,; — 0, 6; and 
05 — +00, so from Eq. (2.83), we have «419 — 2, and Eq. (2.82) becomes 


SA Ola to 2 [7] — 
A/ug 21-e Vi : »142e Vd (2.84) 
cosh 0, + cos 63 


If the damping ratio, & = 0, the right-hand side of Eq. (2.84) is equal to 3, which 
implies that the maximum pseudo acceleration is equal to 3 times ug. For actual struc- 
tures with damping, the normalized maximum pseudo acceleration, A /ug, is always 
between 1 and 3, as demonstrated by Eq. (2.84). This is because when c, /« — 0, 
the ground motion period is much larger than the natural period of the sliding system, 
so there is sufficient time for the oscillation of the relative displacement to decrease 
and stabilize to the static displacement, u,, = —ôug/œ?, prior to the onset of the 
opposite sliding. Once sliding changes direction, the relative displacement initiates 
oscillation around the new static equilibrium displacement, — usr. If the damping 
is sufficiently small, the relative displacement can approach — 3u,, during the first 
oscillation cycle; if the damping is very large, the oscillation of the relative displace- 
ment diminishes rapidly and is only capable of attaining the new static equilibrium 
displacement, — ust. Figure 2.13 shows the steady state response of the normalized 
pseudo acceleration for a = 0.8, $ = 5%, T = 1 sand T, = 10 s (w/w; = 0.045). 

Equation (2.82) can also be written as 


~—Leowt 1 
A/ug=1+2e ve 3 (2.85) 
(e-9 — 1)" + 2e-? (cos 02 + 1) 


If € = 0 [which means 0; = 0 from Eq. (2.31)] and 6; = (2n — 1)z (where n is 
a positive integer), the normalized maximum pseudo acceleration, A /ug, becomes 
infinite from Eq. (2.85). If & Æ 0, a local peak value will be obtained for A / ug when 
0» = (2n — 1)r, i.e., wia = Qn — 1)o, from Eq. (2.31), in which c; is the natural 
frequency of the damped vibration during the sliding phase: 


* Sliding direction changes 


t (s) 


Fig. 2.13 Normalized pseudo acceleration response for o = 0.8, T = 1 s, Tg = 10 s and £ = 596 
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epee 


A/ug=1+2e v" 


T 1 


—e- 


(2.86) 


A 


This result is consistent with the observations from Fig. 2.12b that there are 
multiple resonant frequencies in the area of w,/a@, < 1. 

For £ = 0, Fig. 2.14 displays the normalized pseudo acceleration response for 
€, = wz and @ = 3«,. The ground acceleration, üg, is taken as a, cos (wet), where 
Gg is a large value so that sliding can occur during the whole time of investigation. 
Since «4 = @ = (2n — 1)o, for $ = 0, the relative displacement reaches the peak 
value of a certain half sliding cycle when the direction of sliding changes, which is 
the furthest location from the static equilibrium displacement of the next half sliding 
cycle. From Eq. (2.13), the relative displacement at the end of one half sliding cycle 
for $ = 01s 


u, (t; + T,/2) = cos(wı Tg /2)us(t;) + us; — usi COS(@ T,/2) = —us(t;) + 2us; 
(2.87) 


So 
lus (t; + 7,/2)| — lu (t;)| = 2ug/o? (2.88) 


Thus, during one cycle of the ground motion, 7',, the maximum relative displace- 
ment can increase by 4g /o, and can continue increasing until the ground accelera- 
tion is not large enough anymore to start sliding. For actual structures, the maximum 
relative displacement will be reached after several cycles of ground motion because 
of the damping, as shown in Fig. 2.15 for $ = 5%, œ = 0.8, T = 1 sand T, = 0.45 s 
(vs / a = 1). When w,/@, — 0, 6; — +00 from Eq. (2.31), so Eq. (2.86) tends 
toward Eq. (2.84). This means that when c, /c decreases, the resonance will slowly 
disappear, as shown in Fig. 2.12b. 

When w,/@,; — +00, 0, and 0) — 0 from Eq. (2.31), we have wiqTo — 65/2 
from Eq. (2.81), so Eq. (2.80) becomes 

1—28 (2 ) 
m 
8 [on 


I 2e: 
A/ug=|l-e vi" 
cosh 6; + cos 62 


This is because when c/c — +00, the frequency of ground acceleration is 
considerably larger than the natural frequency of the sliding system, so the oscilla- 
tion of the relative displacement can be hardly stimulated because of the frequently 
changed sliding direction. 


>0 


(2.89) 
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Fig. 2.14 Response of the normalized pseudo acceleration for £ = 0: a w1 = wg; b @ = 3c, 
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Fig. 2.15 Response of the normalized pseudo acceleration for § = 5%, a = 0.8, T = 1 s and T, 
= 0.45 s 
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2.5.3 Solutions for the Sliding Displacement Amplitude 


In Sect. 2.4.2, we have obtained 


ÔT QO fl. 
cos(o,t;) = — + A 
8 8 
dug AD, sin 05 T 


= + (2.90) 


as \ 2/1 — E? cosh; +cos% 2 


By substituting f = t; + 2/@g, into Eq. (2.18) and using Eq. (2.90), the amplitude 
of the sliding displacement can be obtained as follows: 


us, ap = us(ti + z/Gg) — us (ti)| 


1 ; 2ag . ma 
= gone (Ios) + o(2u,o + à,o(/og)) — z sin(wgt;) — P cos(ogt;) 


8 & 
—2ag , 
= —3 Sin(wgti) + 2auyo 
& 
— sinh 6j + 1 — & Jsine 
E 48 sin(wpt;) +20. i g i) i (2.91) 
we w2 cosh 0; + cos 05 


Therefore, the normalized sliding displacement amplitude is 


T sinh 6; + g 1— sP) sin 05 


2 . g 
s.d = |2 ti 2 
is pl (asi) sin(, ) ee aga cosh 6, + cos 62 


(2.92) 


As revealed by Eq. (2.92), unlike the maximum pseudo acceleration, the sliding 
displacement amplitude is dependent on the amplitude of the ground acceleration, 
Gg, for the sliding-sliding case. When a,/ug — +00, cos(wgti) — 0 as revealed 
by Eq. (2.90), so the normalized sliding displacement amplitude, us,ap/ (ag /o».). 
in Eq. (2.92) tends toward to 2, which is consistent with the results presented in 
Sect. 2.5.2. 
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2.6 Conclusions 


This chapter presents the responses of 2DOF sliding base systems subjected to 
harmonic ground motions. The response history of a SB system can exhibit two 
types of phases: the stick phase and the sliding phase. During the sliding phase, the 
vibration of the relative displacement has a higher natural frequency and damping 
ratio compared with the stick phase, which is related to the mass ratio. The equiva- 
lent dynamic force for the vibration of the relative displacement during the sliding 
phase is a step force, which results a static displacement as given in Eq. (2.12). 
The responses of the sliding base system under harmonic ground motions converge 
rapidly to steady periodic responses with the same period as the ground accelera- 
tion. The sliding base system displays three types of motions, namely stick-stick, 
stick-sliding, and sliding-sliding, as the ground acceleration amplitude increases. 

A sliding base structure has multiple resonant frequencies when it slides, unlike 
fixed base structures. For the sliding-sliding case, resonance happens when the period 
of ground motion is odd times of the natural period of vibration of the relative 
displacement that occurs during sliding. As the mass ratio increases, the maximum 
pseudo acceleration generally decreases. However, Resonance can result in obtaining 
larger maximum pseudo acceleration for larger mass ratios in some particular cases. 
Moreover, in certain instances, the sliding base structure can have a higher maximum 
pseudo acceleration than its fixed base counterpart due to extra resonances. In the 
sliding-sliding case, the upper limit of pseudo acceleration is reached, regardless of 
the amplitude of the ground acceleration. Equations (2.80)-(2.83) can be used to 
calculate the maximum pseudo acceleration in the sliding-sliding case. 

The amplitude of the sliding displacement is a suitable response quantity to repre- 
sent the extent of sliding, and it can be normalized by the half of the vibration 
amplitude of the ground displacement, a, / o». The normalized sliding displacement 
amplitude, us ap / (as /o».). is more affected by the mass ratio and frequency ratio, 
@,/@, in region w/w > 1 compared to region c/c < 1. The normalized sliding 
displacement amplitude exhibits a higher value as a/g increases following the 
sliding occurring, yet its rate of increase gradually reduces, and it eventually reaches 
an upper bound of 2. 
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Chapter 3 A) 
Response Histories of Sliding Base m 
Structures Under Earthquake Excitation 


3.1 Equations of Motion 


Figure 3.1a shows a schematic diagram of an N-story SB structure, the sliding base 
of which rests on the foundation. By lumping the structure mass at the corresponding 
floor level, the multistory SB structure can be analyzed using the simplified model 
shown in Fig. 3.1b. The lateral displacements of the ith floor in the x and y directions 
with respect to the sliding base are denoted as u,4; and u,,;, respectively. By further 
assuming that the structure is symmetrical about the x and y axes, the equations of 
dynamic equilibrium of the structure subjected to three-component excitations are 
as follows: 


m| (üx + iisy)1 + ii, | +c,u,-, + k,u,, = 0 
mp (is, + ds.) + 1'm[ (tier + üsx)1 + ün] = f 
m (icy + üsy)1 + try] + cù, +k,u,, = 0 

) [(ii Ugy + iigy)1 +ü] = fy 


(3.1) 
mp(iigy + lisy) + 1m 


where m, c, (c,), and kx (k,) are the mass, damping, and stiffness matrices, respec- 
tively, when the structure base is fixed; m; is the mass of the sliding base; 0 and 1 are 
the vectors whose elements are all zero and unity, respectively. i x(t) and i,, (1) are 
the x and y components, respectively, of the ground acceleration; tsx (t) and i/;, (1) are 
the sliding accelerations with respect to the ground in the x and y directions, respec- 
tively; uj, = [ur Urx2; .-. Ur N] and u,, = [uryi Upy2 0-5 Uryn | are the floor 
displacement vectors in the x and y directions, respectively; à, and U,,, and ü, and 
ü,, are the corresponding velocity and acceleration vectors, respectively; and f(t) 
and f ,(t) are the x and y components, respectively, of the friction force at the sliding 
interface. In Eq. (3.1), the first and second (third and fourth) equations represent the 
dynamic equilibrium of each floor mass and the whole structure, respectively, in the 
x-direction (y-direction). 
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Fig. 3.1 N-story SB structure: a schematic diagram; and b simplified model 


For the model shown in Fig. 3.1b, the mass matrix m is 


mMm 0. 0 

m = mae 0 (3.2) 
sym : 
my 


where m; (i = 1, 2, 3, ..., N) is the mass of the ith floor. The stiffness matrices k, 
and k, are 


Ka F kyo —ky 0 ib 0 0 
kyo F ka C Kx3 0 0 
kx3 + ky4 + 0 0 
k, = EA : : (3.3) 
—kxw-1) 0 
sym. ky p) + kyn ken 
kyy 


and 
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ky c ko —ky2 0 0 0 
kyo T ky3 —ky3 0 0 
ky3 + kyy 0 0 
k, = PEE : G4) 
—kyw-1 0 
sym. kyw-1) + kyn —kyn 
kyn 


where k,; and ky; are the lateral stiffness of the ith story in the x and y directions, 
respectively. The damping of the superstructure can be assumed to be of the Rayleigh 
type; thus, when the damping ratios of two specified modes of the superstructure are 
given, the damping matrix c, and c, can then be determined. 

The response history of an SB structure can exhibit two types of phases: stick 
phase and sliding phase. For the stick phases, during which sliding does not occur, 
the sliding acceleration is equal to 0, and the friction force, f, is smaller than the 
static friction force; thus, we have 


üs, = lisy =0 
F (3.5) 
f= Je + F <(m+ mp)(g + lige) Ms 


where ii. (t) is the z (vertical) component of the ground acceleration, g is the gravity 
acceleration, and ms is the static friction coefficient. Combining Eqs. (3.1) and (3.5) 
leads to 


m[i,.1 F ii, | F Crù,x F k;u,, =0 
m[i,,1 + ü,,] T 6,0, + ku, = 0 


2 2 
1Tmiü,, 7 1Tmii,, . " 

—— —— — ^ lig; +| ——— + ley] < (8 + gc) 

(= m; + Mp B »» m; + Mp i (s i i 


(3.6) 


The first and second equations of Eq. (3.6) govern the response of an SB structure 
during the stick phases. The third equation of Eq. (3.6) is the precondition for the 
stick phases; when it is no longer satisfied, a sliding phase starts. 

During the sliding phases, the direction of the friction force is opposite to the 
direction of the sliding velocity; thus, we have 
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N 


dp D +m) (8 + tig) mM 
V usx t US$, Niel 


" (3.7) 
—Usy T 
fh = m; + mp |(g + ügz)u 
i JL +i, (x ) a 


where u is the dynamic friction coefficient. Combining Eqs. (3.1) and (3.7) yields 


ii,,m1 T mü, F C, F k,u,, = —ii,, m1 
.. 1"mü,, sx .. E 
usx + N + (g + ügz)u = —Ugx 
Dizimi +m, fi, +12, 
T " " (3.8) 
üsyml + mü,, + c,ü,, + k,u,, = —i,,ml 
1'mii,, ls, 


lisy F F - (g s jigs) = —ligy 
Vier ty dul + 02, 


Equation (3.8) governs the response of an SB structure during the sliding phases. 
When the sliding velocity during a sliding phase becomes 0, this round of sliding ends. 
Afterwards, the structure may continue to slide or enter a stick phase depending on 
whether the third equation of Eq. (3.6) is satisfied. Based on the governing equations 
for the stick and sliding phases and the transition conditions between different phases, 
as presented above, a program was developed for computing the responses of SB 
structures subjected to three-dimensional excitations. The numerical methods for 
solving Eq. (3.8) will be presented in the next section. 


3.2 Numerical Computation Methods 


The first and second equations of Eq. (3.6), which govern the response of an SB 
structure during the stick phases, are the same as the governing equations of a fixed 
base structure. The numerical methods for solving these equations can be found in 
Chopra (2001). 
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The system of differential equations presented in Eq. (3.8) can be solved by the 
time-stepping method. The response quantities (is»)j+1, (s) 4 Oba, Urx) 41 


(Ür) (Uy); p (ry), ,, and (ii,y), ,, at time i + 1 satisfy Eq. (3.8) at time i + 1, 
i.e., 

Gi); ml + mM(ü,x);p1 + 60r); + Kr (ui); = — (iex); ml 
1Tm(ü,,);,, E (üsx)i+1 


ym + y 1 Í 
Xu m Tib O F [oka] 


É + (ites) Jo = —(iigz) 44 
(iisy);,,m1 + m(ü,,), ,, + Cy (ùy); tk, (uy) ic) = —(ügy); ml 


1Tm(ü, 
(iue ud 


(üsx)i41 + 


(3.9) 


i (i), 


N m 2 
Dia Mi + Mp (est as OA 


(g + üe:)u = — (is); 


Using Newmark’s equations (Newmark, 1959) for the relationships between the 
response quantities at time i + 1 and the corresponding quantities at time i, we have 


. : lus lu 
USE = (üsx); + g^t sx); F zA sx) 


(uiii = (Ux); us At (Ùx); =P [(0.5 m B)(At)] (ix); + [BAN] üi 


1 1 
(Ùx); = (à,..); + 5 At rx): + zA Ü) 
: : locus l, ja 
(is) = üs), 21005), 2A), 
(2), = (Wo), + Ari), +105 — DADI»), + BAD] ëo) 


: E cons l. du 
(à,,);., = (iy); + 5^! (ü.,), * z At (üs) 
(3.10) 
where At is the time step, and f, which ranges from 1/6 to 1/4, is a parameter that 


controls the variation in the acceleration over a time step. Equation (3.10) can be 
converted to 
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as 27. : s 
Usx)igi = [asx int = (sx); | — Usx); 


T At 
E NM 1 : 
Gina = ET — (u);] 3 A; Č; (; 3B Jo 
1 ie ae AP 
(Ure )i41 = 25A; ein m (u,x); | + (1 m =) (ù); + (ar E aa) iss 
(i). = [Ga 7 Gn): ] - ëo): 


(à). = gapla- Cod.) 7 gartd- (ag 7 1) 693, 
+ 


(à), = za [9 7 Woh] + (17 25) lio), ar- E) G, 


(3.11) 
Substituting Eq. (3.11) into Eq. (3.9) gives 
ES nid m dien 
T Ulsx Ji x x [Mig] = 
Ar iH BAN 26At SU UN 
2 on (sx); is 
A; 3i uc — 5 ls + (is), Jo 
: 2 ; 
MONS F [CM 
1 1'm(u,~) 4 
2 WN = p2 
B(Ar) 20 m; + my 
2 1 1 (3.12) 
EE [zm tgar” t J Qs) = e 
23. (isy); . 
At (isy); 1 + 7 = 7 E * (is), Ju 
TEN T [69] 
1 1Tm(u,,), 
( xd = p4 


BAD? YI mi + m, 


where 
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: Oe ee E . 
Pi = (ügx); 4 t A; sedi + (us); {m1 4 pane DBAt* (Urx)j 


1 1 . 1 At is 
t | m (: Je: eos ! ( 1)m (^ Je fëm: 
BAt 2B 2p 4B 


E PS " 
pa-— —(iige) 44 s Ar xi + sx)i 


jin | l x); + n+ (5 Jr x] 
EN mj +m LBAN I par i NOB Tx 


és 2. ps " 1 1 
p3 zl (sy); 4 3 A (sy); 4 (is), om Frou za Je 


-[ gsm (1 zole» + | (55 1)m (ar Foën), 


9] 2 >. .. 
pa = — (iigy); 4 4 Ar (sy); + (sy); 
1Tm 


1 1, ) (1 : 
emen rs (8); * BAt (à); * tr JEN 


When determining the response quantities at time i 4- 1, the response quantities 
at time i are already known; thus, Eq. (3.12) is a system of nonlinear equations with 
four unknowns, namely, (;,); 4, (Urx)iai> (sy) id and (u,y), " These nonlinear 
equations can be solved using the Newton-Raphson iteration technique (Chopra, 
2001). After determining the values of (2,,);,4, (Urx)j41> (45,); and (u) 


other response quantities, namely, (u;;); 41, Gisc); 41, Qiu. Ux) j41, (usy) 


(3.13) 


i+1 


icr 
(iisy), "HT (ùy); 4, and (ii,y), ,, attime i + 1 can also be determined using Newmark’s 
equations. 


3.3 Response Histories Under Earthquake Excitation 


Figure 3.2 shows the x-component of the ground accelerations recorded at the 
Mammoth Lakes station during the 1980 Mammoth Lakes earthquake. The peak 
ground acceleration (PGA) in the x-direction is 3.8 m/s”. Figure 3.3 shows the 
absolute acceleration history of the top mass of a single-story fixed base structure 
subjected to this ground motion. The natural periods of the superstructure along the 
two principal axes are both 0.3 s, and the damping ratios are taken as 5%. The peak 
absolute acceleration of the top mass is 7.66 m/s”, which is approximately 2 times the 
corresponding PGA. Figure 3.4 shows the responses of SB structures with different 
friction coefficients subjected to this ground motion. The superstructure is the same 
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Fig. 3.2 x-component of the ground accelerations recorded at the Mammoth Lakes station during 
the 1980 Mammoth Lakes earthquake 
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Fig.3.3 Absolute acceleration history of the top mass of a single-story fixed base structure 
subjected to the Mammoth Lakes record 


as that of the aforementioned fixed base structure. The mass of the sliding base is 
3/7 of the top mass. As shown in Fig. 3.4, the peak absolute acceleration of the top 
mass decreases quickly as the friction coefficient decreases; when u = us = 0.05, 
the peak absolute acceleration of the top mass is 1.16 m/s?, which is only 1596 of 
the peak absolute acceleration when the base is fixed. The response history of the 
sliding displacement varies significantly when different friction coefficients are used. 
As the friction coefficient decreases, sliding occurs more frequently. Furthermore, 
the maximum sliding displacements may be obtained in the opposite directions for 
different friction coefficients. 
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Fig. 3.4 Response histories of SB structures with different friction coefficients subjected to the 
Mammoth Lakes record 
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Chapter 4 A) 
Peak Superstructure Responses geai 
of Single-Story Sliding Base Structures 

Under Earthquake Excitation 


4.1 Critical Parameters and Their Ranges 


In Chap. 3, the equations of motion for multistory SB structures have been derived. 
For single-story SB structures (Fig. 4.1), the governing equations can be simplified 
to the following forms: 


" " 2 “ 
Ury + 2€ Wx Uy F Wy Urx = —Ugx 


ii, + 2Ey@yitpy + O Ury = —iigy (4.1) 


V (ais. + üg) + (aii;y xm ügy) < (g F lige) Ms 


for the stick phases, and 


liss + rx + 26, @yUpy + WM Urx = —Ugx 
ee Usx .. .. .. 
Usx + 3 2 (g + lig.) + Uy, = —ügx 
u$. Usy 
(4.2) 
lis + ü,y + 26y@yUyy + Wy Ury = —Ugy 
Usy 


lisy + (g + iig:)M + Ql py = — ügy 


PP 72 
usx 3E Usy 


for the sliding phases. In Eqs. (4.1) and (4.2), o, = /k,/m and éx = c,/(2ma,) 
(@y = Jky/m and $, = cy/(2may,)) are the natural frequency and damping ratio, 
respectively, of the corresponding fixed base (FB) structure in the x-direction (y- 
direction), and œ = m/(m + mp») is defined as the mass ratio. 

As revealed by Eqs. (4.1) and (4.2), the response of an SB structure is greatly 
affected by the friction coefficient; thus, it is meaningful to relate the intensity of 
the ground motion and the acceleration response of the superstructure to the friction 
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Fig. 4.1 Single-story SB structure: a schematic plot; and b simplified model 


coefficient, u. To consider the response in the x-direction, we introduce 


= 2 (4.3) 


which is the displacement of the corresponding FB structure when subjected to a 


static force, mgu, in the x-direction. Dividing Eqs. (4.1) and (4.2) by u;,, leads to 
the following equations: 

E RN 2— 2 agxo 

Uy, + 25,0, Up + QyUyy = —W, Ugx 

Us 
2 
ae Oy V. w = Agx0 Agy0 + 
u 2E,o,| — Ju aja 
ry + 28) (2) ry + (3 ry "ug aa P (4.4) 


2 2 " 
uis Agy0 = u 2H 
) + (si. + DLE) « (1 + we) og 


and 


m .. si a Os 
es ES j= 2 Agx0 + 
Usy + Urx + PONTEM + Wry, = — 0} ee 


= Usy ü z 2 T uid 
Usx + ——[i + JE + xus = — W Ugy 


(4.5) 


2 
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x 2: y 28x &y 
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ze T 2 
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iy + (14E cs en T Uu. = —0, mp Ugy 
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where u,,(f) = u,.(t)/u;,, and u,,(t) = u,,(t)/u$,, are the normalized relative 
displacements and usx (t) = = Usy (1)/Usex and usy (t) = = Usy (u$. are the normalized 
sliding displacements; u, (t), U,y(t), Usx(t) and u,, (f) and u,, (t), i, y(t), as (t) and 
Usy(t) are the corresponding normalized velocities and accelerations, respectively; 
agxo and agyo are the peak values of the x and y components, respectively, of the 
ground acceleration; and ug (f) = = ligx(t)/dgxo and u5, (t) = iis, (t)/agyo represent 
the waveforms of the ground acceleration history. Equations (4.4) and (4.5) imply 
that the normalized displacements, namely, u,, (f), u,, (t), Usx(t) and Usy(t), are 
dependent only on wy, @y/@x, Ẹx, y, a, Bs Is Agx0/ U8, Agy0/Agx0, Ugx(t), Ugy(t) 
and tg-(t)/g, among which dgyo/dgx0, Ugx(t), ug, (f) and iig: (f)/g are associated 
with the ground motion characteristics. 

The maximum earthquake force applied to the superstructure in the x-direction 
can be expressed as mA,, where A, is the peak pseudoacceleration (Chopra, 2001) 
in the x-direction, which is defined as 


A, = o? x max(|u,, (t)]) (4.6) 


By using Eq. (4.3), the normalized peak pseudoacceleration, A,/1g, can be written 
as 


A, E wr x max(|U;x(t)]) E max (|; (t)|) = max(|i,,(t)|) (4.7) 


HE HE Ustx 


Therefore, Ax/ug is equivalent to the maximum absolute value of the 
normalized relative displacement, w,,(r). Since w,,(f) is dependent only on 
Wx, Oy / Oy , Ex, £y, Q, Iis / Lo, agx0/ U8, agy0/4gx0, Ugy (t). Ugy (t) and lis; (t)/g, Ang 
is also only dependent on these parameters. According to the principle of symmetry, 
the parameters that determine the response in the y-direction are the same as those 
that determine the response in the x-direction. 

The first step towards conducting parametric studies in the following sections 
is to investigate the ranges of the critical parameters in accordance with common 
practice. Due to the fact that SB structures are designed for low-rise buildings, the 
natural periods of the superstructure, denoted as Ty = 27/@, and T, = 2zt/oy in the 
x-direction and y-direction, respectively, are limited under 1.0 s. Building structures 
are typically designed with similar stiffnesses in two orthogonal directions; therefore, 
the range of T/T, = w,/a,, the ratio of superstructure periods in the two horizontal 
directions, is assumed to be from 1/2 to 2. For the general application of SB isolation 
in masonry structures (e.g., Nanda et al., 2015; Qamaruddin et al., 1986), bond beams 
are constructed under the masonry walls as the SB element; since the bond beams 
weigh less than the roof (or floor) diaphragm, the resulting mass ratio, a, will be 
larger than 0.5 (Qamaruddin et al., 1986). When using sliding isolation bearings 
(e.g., Jampole et al., 2016), the mass of every floor is considered nearly equal; thus, 
the resulting mass ratio is approximately 0.5 for single-story buildings and over 0.5 
for multistory buildings. Based on the statements provided, the mass ratio, o, is 
taken to be not less than 0.5 for subsequent analyses. Several studies (Barbagallo 
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et al., 2017; Nanda et al., 2012, 2015; Yegian et al., 2004) investigating friction 
characteristics of sliding interfaces have found that the static friction coefficient, us, 
is slightly greater than the dynamic friction coefficient u; previous studies (Yegian 
et al., 2004) found that the maximum observed value of j15/j4 was 1.38. Therefore, us/ 
u is taken to range from 1.0 to 1.4. The damping ratios in the x and y directions (i.e., 
€, and £,, respectively) are both taken as 5%, which is a commonly adopted value. 
The range of the values of the dynamic friction coefficient, u, of the sliding interfaces 
investigated for SB structures (Barbagallo et al., 2017; Hasani, 1996; Jampole et al., 
2016; Nanda et al., 2012; Qamaruddin et al., 1986; Yegian et al., 2004) is from 0.07 to 
0.41. Furthermore, since the peak ground acceleration (PGA) seldom exceeds 1.2g, 
the maximum value considered for agxo/ ug is taken as 20, i.e., agxo = 1.4g if u = 
0.07. 


4.2 Earthquake Ground Motions Considered 


The ground motion records were selected for each site class defined by ASCE 7- 
10 (ASCE, 2010) from the Pacific Earthquake Engineering Research Center-Next 
Generation Attenuation (PEER-NGA) database. Since SB isolation takes effect 
primarily under significant ground accelerations, it is advisable to select ground 
motion records that have sufficiently large PGAs. by doing so, an extremely small 
value of u can also be avoided as far as possible when conducting a parametric study. 
Therefore, for site classes C and D, only the ground motion records in which the peak 
value of the x-component, agxo, is not less than 0.152 were considered. However, less 
than 50 records with agxo > 0.05g could be found in the PEER-NGA database for site 
classes B and E; as a result of this limitation, 40 acceleration records were selected 
for site classes B and E, each with agxo > 0.05g, and all of them were non-pulse-like 
records. Site classes C and D each had 120 acceleration records selected; within 
those 120 records, there were 90 non-pulse-like records and 30 near-fault pulse-like 
records in each group. The 90 non-pulse-like records for both site class C and site 
class D were selected according to different combinations of the magnitude interval 
and source-to-site distance (defined as the closest distance to the fault rupture zone) 
interval. Three intervals of the magnitude, M, namely, 5.2 € M «6.0,6.0 € M «6.7 
and 6.7 < M < 7.7, and three intervals of the source-to-site distance, D, namely, 0 < 
D «14km, 14 < D «24 km and 24 < D < 120 km, were employed herein, resulting 
in 9 different combinations of magnitude and distance intervals. Given a specific 
site class (C or D), 10 records were selected for each of these combinations. A total 
of 320 records of earthquake ground motion, originating from 69 earthquakes with 
a magnitude M ranging from 5.2 to 7.7, were used in this study. Figure 4.2 shows 
the distribution of the magnitudes and source-to-site distances of the ground motion 
records selected for each site class. 

The previous section mentioned that the response of the superstructure in the 
x-direction is influenced by both the y and z (vertical) components of the ground 
acceleration; therefore, comparing the PGAs of all three components is necessary to 
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Fig. 4.2 Distribution of magnitudes and source-to-site distances of the ground motion records: 
a site class B; b site class C; c site class D; and d site class E 


investigate the significance of the interaction between the two horizontal directions 
and the impact of the vertical component. Figure 4.3 shows the distributions of a5,9/ 
agxo ANd Ay-0/Ag,9 for the 320 selected ground motion records. As shown in Fig. 4.3a, 
the values of a,,9/a5,9 are mostly (95.6%) between 0.5 and 2 with an average value of 
1.02. This shows that the assessed earthquake ground motions typically demonstrate 
comparable intensities in both horizontal directions. Figure 4.3b reveals that the ag-o/ 
agxo Values are primarily concentrated within the 0.2—0.8 range, with an average value 
of 0.62, which indicates that the PGA of the vertical component is generally smaller 
than those of the horizontal component. 


4.3 Normalized Ground Motion Intensity for the Initiation 
of Sliding 


If the third equation of Eq. (4.1) is satisfied throughout the entire excitation history, 
sliding will not occur; thus, the critical static friction coefficient, 4er, for the initiation 
of sliding can be determined using the following equation: 
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where ii,,(f) and ij,,(t) are computed using the first and second equations of 
Eq. (4.1). According to Eq. (4.8), disregarding the influence of the vertical 
component, i/;.(f), there exists a linear correlation between pe, and PGA; and 
hence, max (agxo, agyo) /ug can be used as a more generalized indicator for deter- 
mining the occurrence of sliding, i.e., sliding occurs when max (agxo, agy0) [Usg > 
max (agxo, agy0) / bcr 8 

Figure 4.4 plots the counted median, counted 5th percentile, and counted 95th 
percentile values of max (a gx0, 4 2»0) / Ler & for the 120 selected ground motion records 
for site class D. In the computation, T, is taken to be the same as T,. Figure 4.4a 
shows that as the mass ratio, o, decreases, the median value of max (agxo, dgyo) / cr 8 
increases for a given superstructure period, Ty. The reason for this is that, in 
short-period structures, the peak relative acceleration is typically larger than the 
corresponding PGA, leading to larger values of uer in Eq. (4.8) for larger values 
of a. Figure 4.4a also shows that, except for Ty < 0.3 s, the median value of 
max (agxo, a gy0) /Ucrg increases as T, increases for a given mass ratio. This result is 
consistent with the shape of the corresponding response spectrum of FB structures. 
When max (agxo, agy0)/ Msg is equal to the corresponding counted 5th percentile 
value of max (a gx0, 4 20) / Ler 8, Sliding occurs for a small number of ground motions. 
However, the effect of this short-term sliding on the superstructure response is 
insignificant because of the very short sliding duration in these cases. Therefore, the 
response of an SB structure can be considered the same as that of the corresponding 
FB structure when max (agxo, agy0) /usg is smaller than the corresponding value in 
Fig. 4.4b. During 0.2 s < T, < 0.4 s, the 5th percentile value of max (agxo, agyo) /Mer& 
is estimated to be 0.27 for a = 0.9, and it increases to 0.42 when o = 0.5. The values 
in Fig. 4.4c can be regarded as the lower bounds of max (agxo, Agy0) / [Ls g to ensure the 
occurrence of sliding. The trends observed in Fig. 4.4 are also evident in the results 
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Fig. 4.4 Typical statistical values of max(dgx0, dgy0) /Acr8: a counted median; b counted 5th 
percentile; and ¢ counted 95th percentile 


for the other site classes, except with slight difference in their respective specific 
values. 


4.4 Parametric Study for the Maximum Superstructure 
Response 


According to the previous discussion on the range of dg,0/jug, ten levels of agxo/ Ug, 
namely, 0.25, 0.5, 1, 2, 4, 6, 8, 12, 16 and 20, in which dg,o/ug = 0.25 is basically 
equivalent to the FB case, are used in the following analyses. The dynamic friction 
coefficient, u, for the sliding interfaces used in SB structures (Barbagallo et al., 
2017; Hasani, 1996; Jampole et al., 2016; Nanda et al., 2012; Qamaruddin et al., 
1986; Yegian et al., 2004) falls between 0.07 and 0.41. If u is limited to the range 
between 0.07 and 0.41, none of the selected ground motion records can yield all 
the levels of aşxo/ug considered. However, in order to analyze the dispersion of the 
superstructure response at various levels of d,,o/jwg, it is necessary to apply the same 
number of ground motion records for all a,,0/1.g levels. For this purpose, the value of 
u is adjusted with unscaled ground motion records for each target value of a4:0/1.g. 
By doing so, to reach a large value of agxo/ ug (e.g., agxo/ug = 16 or 20) for ground 
motion records with small PGAs, it will be inevitable to use very small values of u 
(e.g., u < 0.02). 

Figure 4.5 shows each individual value of A,/jzg computed using the 90 non- 
pulse-like records for site class D in addition to the mean, the mean plus one standard 
deviation (SD) (corresponding to the 84th percentile value of the normal distribution), 
the counted median and the counted 84th percentile. In this computation, us = u and 
Tx = T, = 0.3 s are adopted. The mean and the counted median agree well with each 
other, as do the mean plus SD and the counted 84th percentile. As stated above, the 
value of u is adjusted to reach the target value of a4,9/ug with unscaled ground motion 
records. The considered cases are classified into two groups based on the values of 
u obtained; one with u within the common range of 0.07—0.41 and the other with u 
outside this range. Different symbols are used in Fig. 4.5 to denote the data in these 
two groups. As expected, the resulting values of u are basically beyond the common 
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Fig. 4.5 Individual and some statistical values of A, /ug: a æ = 0.5; and ba = 0.8 


range when a4,90/4g = 12, because there are very few records of PGA > 0.8g in the 
ground motion database. In order to determine the reliability of the results obtained 
using u which is out of the common range, the probability densities of the computed 
values of A,/ug at agxo/ug = 1 and 4 corresponding to the two different groups are 
compared in Fig. 4.6; in this figure, the normal probability density functions with the 
corresponding mean and SD are also presented. The distributions of A, /ug in each 
group are fundamentally similar to each other. Based on this observation, similar 
results would likely be obtained for large values of a4:0/1.g if a sufficient number 
of ground motion records with large PGAs were used. Additionally, the probability 
density of the calculated values corresponds quite well with the corresponding fitted 
normal probability density function, which suggests that the probability distribution 
of A./ug for a given agx0o/ug value is approximately in accordance with a normal 
distribution. 


4.4.1 Comparison of the Response in Two Orthogonal 
Directions 


Figure 4.7 shows the mean values of A,/j1g at each level of agxo/ug, in addition 
to each individual value of A,/jg at the corresponding level of agyo/ ug computed 
using the 90 non-pulse-like records selected for site class D. In this computation, s 
= mand T, = T, = 0.3 s are used. Figure 4.7 shows a basically uniform distribution 
of the discrete points with (agyo [RE; Ay/ Lg) coordinates along both sides of the 
mean A,/g versus agxo/ug curve. This indicates that the normalized peak pseu- 
doacceleration and the normalized PGA have an essentially identical relationship in 
both orthogonal horizontal directions; in other words, the results obtained for the 
x-direction can also be applied to the y-direction. Therefore, only the response in the 
x-direction is analyzed hereafter. 
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4.4.2 Effect of the Vertical Earthquake Component 


The responses of SB structures under only the two horizontal components of earth- 
quake excitation were also computed in order to study the effect of the vertical 
component. Figure 4.8 shows the ratios of A, under three-component excitation to 
that under the corresponding excitation with the two horizontal components. The 90 
non-pulse-like records selected for site class D with us = u and T, = T, were used 
to compute the results shown in Fig. 4.8. According to these figures, the vertical 
component of ground motion can either increase or decrease the horizontal response 
of the superstructure. In general, the vertical component has a greater effect on stiffer 
structures. For certain ground motions, with T, < 0.3 s, the ratios of A, under three- 
component excitation to that under two-component excitation exceed 1.15. However, 
for most of the ground motions considered, the ratios of A, under three-component 
excitation to that under two-component excitation are between 0.95 and 1.05, and 
the mean values are basically equal to 1.0. Therefore, the overall effect of the vertical 
component on the superstructure response is negligible. 
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Fig.4.8 Ratio of A, under three-component excitation to that under the corresponding two- 
horizontal-component excitation: a a = 0.5, T, = Ty = 0.3 s; ba = 0.5, dgxo/ug = 4; ca = 
0.8, T, = Ty = 0.3 s; and d a = 0.8, agxo/ ug = 4 
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4.4.3 Effect of the Natural Period of the Superstructure 


Figure 4.9 shows the relationship between the mean A,/ug and T, for different 
values of o and dg.9/ug. The mean values of A,/pg were computed using the 90 
non-pulse-like records for site class D while assuming that us = u and T, = Ty. 
According to Fig. 4.9, as T, increases, the mean values of A, /ug increase for T, < 
0.3 s, but decrease for Ty > 0.4 s, resulting in the maximum mean values of Ax /ug 
are obtained at T, = 0.3 s or 0.4 s. In general, the mean values of A,/ug at T, = 
0.3 s or 0.4 s are close to each other. As dg,0/j1g increases, the influence of T, on 
the superstructure response decreases; for example, the ratio of the minimum to the 
maximum mean values of A, /jwg in Fig. 4.9a is 0.69 for dg.9/4g = 2 and increases 
to 0.87 for agxo/ug = 12. As shown in Fig. 4.9, it can be inferred that the mean 
value of A, /jg does not exhibit a considerable variation within the range of T, that 
has been considered; therefore, it is appropriate to use the response from the period 
with the maximum mean A, /j1g as a representation of the responses for possible SB 
structures. As mentioned above, for site class D, this period can be taken as 0.3 s; 
for site classes B, C, and E, the critical periods obtained are 0.2, 0.2, and 0.4 s, 
respectively. 

In all the above analyses, Tą = T, is adopted. Figure 4.10 compares the mean 
values of A,/jg for different values of T/T, in order to investigate the possible 
impact of T ,/T , on the superstructure response. In general, the mean values of A, / ug 
at a given level of a¿+0/ug decreases as T/T, increases. Nevertheless, this variation 
is quite limited; the presented results in Fig. 4.10 show that the ratios of the mean 
value of A,/ug for T,/T, = 0.5 to T/T, = 2 do not exceed 1.06. Therefore, the 
results obtained for Ty = T, are representative of those obtained for the possible 
value of T,/T, in the range considered. 
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Fig. 4.9 Relationship between mean A,/ug and Ty: a œ = 0.5; and ba = 0.8 
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Fig. 4.10 Effect of T,/T, on the mean values of A,/uwg:aa = 0.5; and ba = 0.8 


4.4.4 Effect of the Difference Between the Static 
and Dynamic Friction Coefficients 


Based on the investigation mentioned above, j1s/j generally within the range of 
1.0-1.4. Figure 4.11 compares the mean values of A,/ug corresponding to different 
values of us/u in order to investigate the influence of jz;/u. As expected, for most 
of the ground motions considered, the superstructure response is not significantly 
influenced by the value of 2/2 when dg,o/4g = 0.25, because sliding does not occur 
at this a¿xo/ug level. The influence of 1,,/j is most clearly observed for dgo/ug = 1. 
If agxo/ug exceeds 1, the influence of 14/1, decreases as agxo/ug increases because 
for larger ag¿xo/ ug values, the responses of SB structures are primarily dominated by 
the sliding phase, during which the responses are independent of the static friction 
coefficient. The ratios of the mean value of A,/uwg for uslu = 1.4 to that for u,/u = 
1 are all below 1.08 and mostly below 1.02. Therefore, the effect of the difference 
between the static and dynamic friction coefficients can be neglected. In the following 


analyses, 4s = u is assumed. 


4.4.5 Effects of the Earthquake Magnitude 
and Source-to-Site Distance 


The mean values of A,/ug for each distance interval and magnitude were computed 
using 90 non-pulse-like records selected for site class D in order to study the effects 
of the earthquake magnitude and source-to-site distance. Figure 4.12 shows the mean 
values of A, /jg for the three magnitude intervals. For the entire range of agx0/ U8, 
no general trend can be observed for A, /jg as the earthquake magnitude increases. 
However, when ag¿x0/ug = 0.25, the mean values of A,/jg for 6.7 < M « 7.7 are 
approximately 10% larger than those for 6.0 « M « 6.7; the relative differences in 
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Fig. 4.11 Effect of j1;/ on the mean values of A,/wg: a æ = 0.5; and ba = 0.8 
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the mean values of A, /ug between any two of these groups at a given agxo/ ug are 
all below 7% and mostly below 5%. This indicates that the earthquake magnitude 


has little effect on the superstructure responses of SB structures. 
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Fig. 4.13 Mean values of A, /jg for three distance intervals: a œ = 0.5; and b a = 0.8 


Figure 4.13 shows the mean values of A,/g for three distance intervals. The 
influence of the source-to-site distance on the superstructure response, like that of 
the magnitude of the earthquake, is also insignificant. 


4.4.6 Effect of Near-Fault Pulses 


The acceleration, velocity, and displacement histories of near-fault ground motions 
influenced by forward directivity contain distinct pulses. To investigate the possible 
effects of these pulses, the mean values of A,/ug computed using the 30 pulse- 
like records and 30 non-pulse-like records both recorded for 0 « D « 14 km and 
site class D are compared in Fig. 4.14. When sliding basically does not occur, 
Le., dgxo/Ug = 0.25, the mean value of A,/jg corresponding to the non-pulse- 
like records is 1.16 times that corresponding to the pulse-like records. The corre- 
sponding FB structure has a larger response amplification factor, A,/a,,9, for non- 
pulse-like ground motions compared to pulse-like ones, as indicated. Chopra and 
Chintanapakdee (2001) reported similar results, where they examined the normal- 
ized response spectra of harmonic excitations containing different numbers of cycles 
to interpret this phenomenon. The response amplification factor increased as the 
number of cycles increased, implying that the response amplification factors of 
pulse-like ground motions with one or several dominant pulses are generally smaller 
than those of non-pulse-like ground motions with more excitation cycles. As illus- 
trated in Fig. 4.14a, when œ = 0.5, the ratio of the mean value of A,/pug for the 
non-pulse-like records to that of the pulse-like records remains almost the same as 
agx0/ ug increases. However, when o increases to 0.8, the ratio decreases signifi- 
cantly as agx0/ ug increases, as shown in Fig. 4.14b. This phenomenon is consistent 
with the fact that the difference in the superstructure response between different 
ground motions decreases as o increases. More detailed results related to this fact 
are presented in the next section. 
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Fig. 4.14 Effect of near-fault pulses on the mean values of A,/ug: aa = 0.5; and ba = 0.8 


4.4.7 Statistical Results for Each Site Class 


Figure 4.15 shows the mean values of A,/jg for each site class. Herein, only the 
non-pulse-like records were used for the computation to maintain consistency. Thus, 
the numbers of ground motion records used for site classes B, C, D and E are 40, 90, 
90 and 40, respectively. As shown in Fig. 4.15, except for agxo/ug = 0.25 (which is 
basically equivalent to the FB case), as the mass ratio, œ, increases, the superstructure 
response reduces. This phenomenon can be easily explained by using the governing 
equations for the sliding phases during unidirectional excitation. Under unidirectional 
excitation in the x-direction, Eq. (4.2) can be simplified into 


2 B 
& Cx wy Us pg 


Úrs + —Urx 
Vl-aJ/l—-a l-a 


which is the differential equation of a single-degree-of-freedom (SDOF) system with 
a natural frequency of c /4/1 — o and a damping ratio of $,/4/1 — o subjected to 
a step force corresponding to a static displacement of ug /œ?. The damping ratio of 
this equivalent system increases as o increases, leading to a general decrease in the 
response of u,, as œ increases. 

After dg.o/ug exceeds 0.5, the mean A,/pg versus agxo/ug curves experience 
a rapid decline in tangent slopes as a result of sliding. When dg,9/jg exceeds a 
sufficiently large value, i.e., there is an upper limit for the superstructure response of 
an SB structure, the tangent slopes are expected to finally become 0. This situation is 
favorable for isolating extremely large earthquakes. The efficiency of the SB system 
can also be demonstrated by the value of A,/a,,9, which is equal to the origin-oriented 
secant slope of the A,/j1g versus agx0/ug curve. Taking œ = 0.8 in Fig. 4.15b (site 
class C) as an example: when aṣx0/ ug = 2, the mean value of A,/a5,9 is 0.82, whereas 
this value is 2.23 for the FB case; consequently, the superstructure response of the 
SB structure is just 36.8% of that of the corresponding FB structure in this instance. 


ü +2 (4.9) 


|usx| 1— æ 
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Fig. 4.15 Mean values of A,/g for each site class: a site class B; b site class C; c site class D; 
and d site class E 


Figure 4.16 presents the ratios of the mean value of A,/jg for site class C (D or 
E) to that for site class B, for o = 0.5 and 0.8, in order to investigate the effects of 
local site conditions on the superstructure response. It is evident that the response 
of the superstructure is affected by the local site conditions; as the site soil becomes 
softer, the mean values of A,/j1g increase. For agxo/ ug = 0.25, when sliding basically 
does not occur, the ratios of the mean values of A,/jwg for site classes C, D and E 
to that for site class B are equal to 1.19, 1.20 and 1.24, respectively. These ratios 
generally decrease as agxo/ ug increases. As shown in Fig. 4.16a, when dg.0/ug = 2, 
these ratios decrease to 1.07, 1.07 and 1.12 for site classes C, D and E, respectively; 
and they further decrease to 1.02, 1.03 and 1.04 when a4,,90/ug = 12. 

Figure 4.17 shows the coefficients of variation (COVs) of A,/jg for every site 
class in order to investigate the dispersion of the superstructure response at a specified 
value of dgx9/ug. The COVs of the A,/wg versus agxo/ ug curves are similar for all 
four site classes. In the range of agxo/ug < 2, the COVs of A,/ug decrease rapidly 
as agxo/ ug increases; after agxo/ ug exceeds 2, the COVs of A,/ug are quite constant. 
This means that sliding tends to reduce the dispersion of the superstructure response 
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Fig. 4.17 Coefficients of variation of A, /jg for each site class: a site class B; b site class C; c site 


class D; and d site class E 


due to the record-to-record variability. As œ increases, the COVs of A, /ug decrease 
in general. This can also be interpreted by using Eq. (4.9); larger values of a leads to 
larger equivalent damping ratios for the sliding phases, further resulting in smaller 


dispersion of the structural response. 
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4.5 Simplified Equations for Estimating the Maximum 
Superstructure Response 


For the design of SB structures, it is desirable to employ simplified equations to 
estimate the peak pseudoacceleration of the superstructure. Based on the preceding 
discussions, we know the following: (1) the relationship between A, / 11g and agx0/ Mg 
is basically identical to that between A,/4 and a,,9/14g; (2) the dependencies of 
the mean value of A, /jg (or A,/jg) on the vertical earthquake component, T/T), 
Uslu, the earthquake magnitude and the source-to-site distance can be neglected; 
and (3) the response of possible SB structures can be represented by the response 
at a critical period for each site class with appropriate conservativeness. Thus, the 
following equation was proposed to estimate the mean values of A, /ug and Ay/pg: 


Ax — Bi (aexo/ ug)” and Ay e By (asyo/ng)^ (4 10) 


H8 — (agxo/ Hg)" + Bs H8 — (ago/ug)^ + fs 


where 61, 62 and £5 are the regression coefficients that depend on the site class and the 
mass ratio, o. Equation (4.10) captures the trend of A, /ug (Ay / ug) with respect to 
agxo/ US (deyo/ Mg); i.e., Ax/u8g — 0 when agṣx0/ ug — 0, and Ax/ug approaches an 
upper limit when d,,9/ug — +00. The Curve Fitting Toolbox of MATLAB (2014) 
was used to conduct nonlinear regression analyses for determining the values of 
the regression coefficients in Eq. (4.10). In order to assess conservativeness, regres- 
sion analyses were conducted using results obtained from non-pulse-like records, 
considering that responses under pulse-like ground motions generally exhibit smaller 
values compared to those under non-pulse-like ground motions. Table 4.1 presents 
the values of 6,, 85 and £4 obtained for each site class and various values of a. The 
values predicted through Eq. (4.10) are also compared to the values computed from 
response history analyses shown in Fig. 4.15. The mean values of Ay/ug can be 
accurately estimated by the proposed equation. The values of 61, 62 and 3 for the 
a values not included in Table 4.1 can be found by linear interpolation of the values 
from Table 4.1, because each regression coefficient has an approximately linear rela- 
tionship with a; if a > 0.9, using the £j, 2 and 63 values for a = 0.9 can give 
conservative results. 

As shown in the previous section, anormal distribution is appropriate for modeling 
the probability distribution of A,/jzg with respect to a specific site class, agxo/ ug and 
a. Therefore, if we can further derive a simplified equation for the COV (or SD = 
mean x COV), then the value of A,/jwg (or Ay/j1g) corresponding to any probability of 
nonexceedance can be readily determined. Figure 4.17 shows that, despite possible 
variations of the exact values among different groups, the trends of the COVs of 
Ax/ug with respect to agxo/ug remain consistent for the four site classes. Because 
the computed COV values are related to the selected ground motion records used 
in the computation (i.e., a different set of records for the same site class may lead 
to different COVs), it is reasonable to expect equivalent dispersion levels for the 
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Table 4.1 Values of the regression coefficients in Eq. (4.10) 


Site class a Bi B2 B3 
B 0.5 3.26 1.07 1.12 
0.6 3.20 0.92 1.22 
0.7 3.10 0.81 1.30 
0.8 2.80 0.73 1.24 
0.9 2.41 0.58 1.15 
C 0.5 3.32 1.03 0.90 
0.6 3.17 0.88 0.97 
0.7 2.92 0.76 0.97 
0.8 2.71 0.63 1.02 
0.9 2.48 0.45 1.10 
D 0.5 3.37 1.01 0.92 
0.6 3.09 0.90 0.90 
0.7 2.93 0.76 0.97 
0.8 2.78 0.61 1.06 
0.9 2.39 0.47 0.99 
E 0.5 3.27 1.16 0.81 
0.6 3.00 1.06 0.80 
0.7 2.73 0.94 0.81 
0.8 2.68 0.68 0.99 
0.9 2.87 0.40 1.41 


four site classes, provided that a sufficient number of records are selected for each 
group. In general, the computed COV values are the largest for site class C under 
the circumstance of the ground motion records considered. Therefore, the data of 
site class C are used to derive the simplified equation for the COVs of A,/ug and 
Ay/ig since a larger COV value leads to a conservative result for a probability of 
nonexceedance larger than 50%. The proposed equation is given by 


5A, /ug = Y1 eXP(—Y2(agx0/Hg)) + ys. and 84, jug = yi €XP(—Y2(aeyo/Mg)) + ys 
(4.11) 


where 64, /yg (B4 jue) is the COV of A./ug (A,/u&). and yı, y2 and y; are the 
regression coefficients that depend on a. Table 4.2 presents the values of y1, y2 and y3 
based on the nonlinear regression analyses for various o values. Figure 4.18 shows a 
comparison of the values predicted using Eq. (4.11) and the values obtained through 
response history analyses. A satisfactory estimation for the COVs of A, / ug is yielded 
by the proposed equation. Similar to Eq. (4.10), for the values of o not listed in 
Table 4.2, the values of y1, y2 and y; can be calculated from the linear interpolation 
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Table 4.2 Values of the 

regression coefficients in E n y2 ¥3 

Eq. (4.11) 0.5 0.23 1.82 0.20 
0.6 0.28 2.25 0.19 
0.7 0.34 2.25 0.15 
0.8 0.42 2.62 0.12 
0.9 0.52 3.23 0.10 


Fig. 4.18 Comparison 

between the COVs of A, /ug 
computed using Eq. (4.11) : 0.5 
0.6 
0.7 
0.8 


and those computed from 
response history analyses 


COV of A/ug 


gy, Hg 


of those provided in Table 4.2; for æ > 0.9, the values of yı, y; and y; for a = 0.9 
can be used. 


4.6 Conclusions 


In this chapter, a comprehensive parametric investigation of the normalized peak 
pseudoacceleration of single-story SB structures subjected to three-component earth- 
quake excitation is presented. The relationship between the normalized peak pseu- 
doacceleration and the normalized PGA is basically identical for the two orthogonal 
horizontal directions. The horizontal response of the superstructure can be either 
reduced or increased by the vertical component of ground motion. If T, < 0.3 s, for 
certain ground motions, the superstructure response can increase by more than 1.15 
due to the vertical component; but the effect of the vertical component is negligible 
for the majority of situations. The normalized peak pseudoacceleration exhibits a 
pattern of initially increasing and subsequently decreasing as the natural period of 
the superstructure increases. For the range of T, considered, the variation in the 
mean value of A,/ug is not very significant. For simplicity and conservativeness, the 
response of possible SB structures can be represented by the response at the period 
where the maximum mean A,/ug is generally obtained. The influence of the natural 
period ratio in the two orthogonal horizontal directions and the possible difference 
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between the static and dynamic friction coefficients on the superstructure responses 
of SB structures is insignificant. 

The effects of the earthquake magnitude and the source-to-site distance are very 
small and can be neglected in practice. Superstructures typically exhibit smaller 
responses when subjected to pulse-like ground motions compared to non-pulse-like 
ones. Local site conditions have an effect on the response of the superstructure. For 
sites located on softer soil, a larger response is obtained, and the dependence on the 
local site conditions decreases as the normalized PGA increases. The trend of the 
COVs of A,/wg with respect to agxo/ ug is similar among the four site classes. The 
COVs decline rapidly at smaller values of a4:0/1.g and remain basically constant after 
agxo/ Ug = 2. For a given site class, a,,0/1,g and g, a normal distribution is appropriate 
for modeling the probability distribution of A,/jg. 

The mean values and COVs of A,/ug decrease as o increases. An upper limit for 
the superstructure response exists for every mass ratio, which is beneficial for the 
isolation of extremely large earthquakes. Implementing Eqs. (4.10) and (4.11) with 
the associated values of the regression coefficients can provide good estimates for 
the mean values and COVs, respectively, of A,/ug and can be used to predict the 
value of the normalized peak pseudoacceleration corresponding to any probability 
of nonexceedance. 
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Chapter 5 (R) 
Equivalent Lateral Forces for Design ES 
of Multistory Sliding Base Structures 


5.1 Model Descriptions 


In Chap. 4, the parameter mass ratio, œ, was introduced. This parameter is defined as 
the ratio between the superstructure mass and the overall mass of the SB structure. 
The response of the superstructure is heavily influenced by the parameter o, as 
demonstrated in Chap. 4. The mass ratio of a multistory SB structure presented in 
Fig. 5.1 can be calculated by 


N 
a= Oinam (5.1) 


La mi + mp 


in which m; is the mass of the ith floor; m, represents the mass of the sliding base; 
and N corresponds to the story number. 

Based on practical applications, a maximum of five stories can be considered as 
the upper limit for the story number of SB structures. For a building with N <5, it is 
reasonable to consider that the mass of every floor is equal. Thus, Eq. (5.1) becomes 


Nm 


= —__ 5.2 
i Nm+my (3) 


in which m represents the mass of each floor. In practical applications, the value of 
a should be at least 0.5 as m is usually equal to or greater than mp. Hence, the value 
of œ can be considered not less than 0.5. 

It is assumed that the stiffness of each story is equal. The parametric study in 
Chap. 4 demonstrates that the natural period ratio between the two orthogonal hori- 
zontal directions has a negligible impact on the superstructure responses. There- 
fore, it is assumed that the story stiffness in both the x and y directions is identical 
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Fig. 5.1 Schematic plot of a multistory SB structure 


and represented as k. Given these assumptions, the fundamental period, 7';, of the 
superstructure is expressed as follows: 


m 
T, =2n fan (5.3) 


in which the values of coefficient C are equal to 0.382, 0.198, 0.121, and 0.081 for N 
— 2, 3, 4, and 5, respectively. As Chap. 4 demonstrates, the superstructure responses 
in single-story SB structures remain basically unchanged when the superstructure 
period falls within the typical range, and it is feasible to conservatively represent 
the responses of potential SB structures with the response at the period that gener- 
ally yields the highest superstructure responses. The truth of this result persists in 
multistory SB structures. The computed results based on T, — 0.3 s are used in the 
subsequent analyses as maximum superstructure responses occur mostly at T; close 
to 0.3 s. 

The construction of the damping matrix involves the utilization of Rayleigh 
damping. For N — 2, the damping ratios for the first and second modes are taken as 
5%, while for N = 3 and 4, the damping ratios for the first and third modes are 5%, 
and for N = 5, the damping ratios for the first and fourth modes are taken as 5%. 

The static and dynamic friction coefficients are assumed to be the same as 
the difference between them was found to have little effect on the superstructure 
responses. In accord with Chap. 4, the maximum value of a,,9/jug is taken as 20, 
where agxo is the peak value of the x component of the ground acceleration, ju is the 
friction coefficient, and g is the gravity acceleration. 
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It was found that the variance between the static and dynamic friction coefficients 
had minimal impact on the superstructure responses, thus they are assumed to be 
equal. In accordance with Chap. 4, an upper limit of 20 is placed on the value of 


Agro / U8. 


5.2 Peak Base Shear 


After obtaining the displacement history by performing response history analysis, 
the equivalent static forces in the x-direction, F, = [Fx1, Fi», ..., Fen] (Fx is the 
force acting on the ith floor) for a multistory SB structure can be determined using 
the following equation (Chopra, 2001) 


F,(t) = kyu (t) (5.4) 
The peak base shear, V,,, can subsequently be computed by 


Vix = max|1"k,u,.(1)| (5.5) 


Because of the close relationship between the peak base shear and the mass of 
the superstructure, as well as the friction coefficient, the normalized peak base shear, 
V bx, is introduced as follows 


u Vs max [17 kyu (£)| 


Vix = —— = 5.6 
"^ Nmgu Nmgu ee 


For a single-story SB structure, Eq. (5.6) is simplified to 


y max|ku,, (t)] o? x max(|u; (t)|) Ax 
bx = i H 
j mgu HE s 


(5.7) 


As indicated by Eq. (5.7), the normalized peak base shear, Vp» is equivalent 
to the normalized peak pseudoacceleration, A, /g, for single-story SB structures. 
The normalized peak pseudoacceleration has been thoroughly examined in Chap. 4. 
Hence, with certain modifications, the equations created for predicting A, /ug could 
be employed to predict Vp, for multistory SB structures. 

By replacing A,/jg in Eqs. (4.10) and (4.11) with V;,, the following equations 
are yielded: 


LA Bi (agr0/mg)” 


Vox (5.8) 
(ags0/mg)” + b3 
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and 
áp, = yı EXP(—Y2(aex0/Hg)) + 73 (5.9) 


in which y,. represents the coefficient of variation (COVs) of V px. In order to verify 
if Eqs. (5.8) and (5.9) are applicable to multistory SB structures, the mean values and 
COVs of V», of structures with varying story number N are compared in Fig. 5.2a, 
b, respectively, using the 180 non-pulse-like ground motion records for site classes 
C and D. As can be seen in Fig. 5.2a, when a,,0/u is given, the mean value of 
V, decreases and ultimately approaches a constant value as N increases. In order 
to account for this impact, the introduction of a reduction factor, denoted as yy, is 
necessary. This factor is calculated by taking the ratio of the mean V, of a structure 
that has N stories, to the mean V», of a single-story structure that has the same 
a under the same a5;0/;4g. The analyzed data suggest that the influence of o and 
agx0/ 4g On yy is minimal. The nonlinear regression analyses yielded the following 
formula for calculating yy: 


yn = 0.25e 9 5N + 0.86 (5.10) 


Therefore, the mean value of V, for an N. -story SB structure can be calculated 
by the following equation: 


f: 
- Agx 
Vix = (0.25e 0" + T ee (5.11) 
(agx0/M8) + Bs 
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Fig. 5.2 Statistical values of Vp: a mean values; b COVs 
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The values of the coefficients £1, 9; and £3 are provided in Table 4.1. As depicted 
in Fig. 5.2b, The COVs of V, are not significantly affected by the story number N. 
Hence, Eq. (5.9) can be used for multistory SB structures without requirement of 
modifications. The values of yı, y2 and y3 are provided in Table 4.2. 


5.3 Equivalent Lateral Force Distribution 


The response histories of V,;/(Nmg) for a three-story SB structure (with a = 0.8 
and u = 0.1) and corresponding FB structure under the Mammoth Lakes record 
(agxo = 0.39g) from the 1980 Mammoth Lakes earthquake are presented in Fig. 5.3, 
where V; is the story shear of the ith story. In the case of the FB structure, the peak 
story shears for various stories occur simultaneously. However, for the SB structure, 
the peak shear time for each story varies. Hence, it is not feasible to employ the 
distribution of equivalent static forces at the peak base shear to ascertain the peak 
shear of other stories. 

Peak story shears are the response quantities that need to be used in design. There- 
fore, the equivalent lateral forces were calculated through the following process: (1) 


Va /(Nmg) 


(a) 


t (s) 


(b) 


Fig. 5.3 Response histories of each story shear of a three-story SB structure (with a = 0.8 and 
u = 0.1) and b three-story FB structure subjected to the Mammoth Lakes record from the 1980 


Mammoth Lakes earthquake 
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Fig. 5.4 Distributions of Ff,/Fy under different percentiles (a = 0.8): a dgxo/ug = 6; 
b dgxo/mg = 12 


Identify the peak story shear of each story for every ground motion; (2) Organize the 
peak story shears for each story from the 180 ground motion records in ascending 
order; (3) Determine equivalent lateral forces using the peak story shears at the corre- 
sponding percentile. The equivalent lateral force of the ith floor is represented as F7. 
The distributions of F7;/ Ffy for various percentiles are illustrated in Fig. 5.4 when 
a = 0.8. As shown in Fig. 5.4, the selected percentile does not have a significant 
impact on the distribution of Ff, / Ffy. Accordingly, the distribution of F7 / Fz, for 
the 50th percentile is selected for the subsequent analyses. 


5.3.1 Parametric Study 


Figure 5.5 depicts the data points (Fį;/ Ffy, i/N) corresponding to different N in the 
same plot. As shown in this figure, the trend is almost the same for the data points 
(Ff, / Few, IIN) for different N. Therefore, the same relationship between F7 / Fry 
and i/N can be used for different N. 

The distributions of F¢,/ Ffy for various a and agxo/ug are shown Figs. 5.6 and 
5.7. The value of F?,/ Ffy for i < N decreases with an increase in œ or agx0/ ug; and 
aS agx0/u8g increases, the distribution of F¢,/F fy becomes fixed. When agx0/ Ug 
exceeds a specific value, the corresponding lateral forces at the lower levels invert 
direction for o > 0.7. These results means that the increase of o or a5,0/,4g results 
in the concentration of equivalent lateral forces at the upper floors. 
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Fig. 5.5 Comparison of the relationships between F°,/F°y and i/N for different N: aa = 0.5; 
bo = 0.8 
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5.3.2 Simplified Equations 


Referring to the shapes of the i / N — F¢,/ Ffy curves depicted in Figs. 5.6 and 5.7, the 
distribution of equivalent lateral forces was modeled using the following equation: 


Fé,/Féy = cG/NY + (0 —0)G/N) (5.12) 


The regression coefficient c in Eq. (5.12) is dependent on a,,0/ug and a. Equa- 
tion (5.12) satisfies the boundary conditions that when i = 0, F7;/ Ffy = 0, and when 
i=N, Fc / Fey = 1. The values of c attained from nonlinear regression analyses for 
varying values of dg,0/j4g and o are presented in Fig. 5.8. The value of c increases 
and approaches a constant as a,,0/44g increases for a given a, and it also increases 
with an increase in o. These trends are in agreement with the effects of aṣxo/ ug and 
o on the distributions of F?,/ Ffy, which has been discussed in the previous section. 
Through nonlinear regression analyses, the following equation was developed for 
calculating the value of c: 


c = (— 4.3a + 0.29)e 9685/48) 4 3.49 — 0.95 (5.13) 
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Fig. 5.8 Values of c for different values of ag.9/ug and a 


As depicted in Figs. 5.6 and 5.7, the distributions of F?,/ Ffy determined using 
Eqs. (5.12) and (5.13) demonstrate good agreement with those that were computed 
based on response history analyses. 


5.4 Conclusions 


In this chapter, the equivalent lateral forces for the design of multistory SB structures 
are studied. The increase in story number N tends to reduce the mean value of the 
normalized peak base shear but has little effect on the coefficient of variation. The 
peak story shears of a multistory SB structure under an earthquake ground motion 
occur at varying times for each story. Based on this truth, the equivalent lateral forces 
can be calculated by utilizing the peak story shears that correspond to the same 
percentile. Based on the computed results, it can be inferred that the distribution 
of equivalent lateral forces is generally not influenced by the percentile chosen. 
The distributions of equivalent lateral forces for various N follow a similar pattern. 
Additionally, with an increase in the normalized PGA and mass ratio, the equivalent 
lateral forces tend to concentrate at the upper floors. Using Eqs. (5.9), (5.11), (5.12), 
and (5.13), the equivalent lateral forces required for the design of multistory SB 
structures can be determined. 
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Chapter 6 A) 
Peak Sliding Displacements of Sliding coset 
Base Structures Under Earthquake 

Excitation 


6.1 Selection of Ground Motion Intensity Measure 


The 180 ordinary ground motion records and 60 near-fault pulse-like records on site 
classes C and D, which has been presented in Chap. 4, will continue to be used in the 
response history analyses conducted in this chapter. In Chap. 3, some typical response 
histories of the sliding displacements are presented. The PSDs in the two principal 
directions (i.e., x and y directions), Us,9 and Usy9, can be directly determined once the 
response histories of the sliding displacements are obtained through response history 
analyses. The maximum of PSDs over all the horizontal directions or the PSD with 
respect to the origin, usro, can also be determined by 


Usro = max y usx (t)? + usy (t)? (6.1) 
t 


In design, the PSDs of interest that are required to check the sliding displacement 
are reliant on the boundary shape of the sliding surface. For SB masonry struc- 
tures (e.g., Nanda et al., 2015; Qamaruddin et al., 1986a), bond beams are generally 
constructed under the masonry walls as sliding elements, and the boundary shape 
of the sliding surface is generally rectangular to align with the building plane. For 
this case, it is necessary to check the sliding displacements separately in the two 
principal directions, and the PSDs required in design are usxo and Uso. Generally, 
when using sliding isolation bearings (e.g., Jampole et al., 2016), the boundary shape 
of the sliding surface is circular. For this case, design requires comparison of usro 
with the sliding displacement threshold. The analyses that follow investigate all of 
Usx0, Usyo, and Uso, covering the two cases mentioned before. 

To develop dependable earthquake excitation-based prediction models for the 
PSDs of SB structures, it is necessary to choose a suitable ground motion intensity 
measure (IM) that provides relatively small record-to-record variability of the PSDs 
at a given IM. The PGA and PGV, which rely exclusively on the ground motion char- 
acteristics, are the most traditional measures of ground motion intensity. Although 
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there have been proven several IMs that consider both ground motion and struc- 
tural properties [e.g., spectral acceleration at the first-mode period of the structure 
(Housner, 1941) and average spectral acceleration (Eads et al., 2015)] to be more 
efficient for seismic response assessment of fixed base (FB) structures, the changing 
dynamic property of a SB structure may render them unsuitable when sliding occurs. 
Thus, a feasible approach to assess the PSDs of SB structures involves choosing an 
IM that solely relies on the ground motion characteristics and analyzing the effect 
of distinct structural properties in isolation. 

PGA and PGV are potential IMs due to their widespread usage among researchers 
and engineers, as well as the availability of corresponding attenuation relationships 
(Villaverde, 2009). To compare the efficiency (Luco & Cornell, 2007) of different 
IMs (i.e., their capability to produce small variability of the PSD at a given IM), 
response history analyses of a SB structure with T, = T, = 0.4 s (T, = 27/a, and 
Ty = 2x[oy), Ex = £, = 5%, and a = 0.7 subjected to the three components of the 180 
ordinary ground motion records were conducted. Figure 6.1 displays the calculated 
values of usxo and uso in relation to their respective PGAs (agxo and agro) and PGVs 
(Vexo and v4:9), where two levels of friction coefficient u are considered, namely u 
= 0.1 and 0.2. In order to maintain consistency with the definition of usro, Eqs. (6.2) 
and (6.3) are respectively used to compute the PGA and PGV corresponding to tso, 
which are the maximum values of PGA and PGV over all the horizontal directions. 


agio = max y iig (t)? + ügy (D)? (6.2) 
Vgro = MAX y tex (£)? + úgy (1)? (6.3) 


The superiority of PGV over PGA is evident from Fig. 6.1. Previous researchers 
(e.g., Jampole et al., 2020; Ryan & Chopra, 2004) have also observed this result. To 
further quantify the performance of different IMs, the results shown in Fig. 6.1 were 
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Fig. 6.1 Correlations between PSDs and corresponding PGAs and PGVs: a u = 0.1; and b u = 
0.2 
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also subjected to correlation coefficient computation; since a linear relationship may 
not be the best representation of the connection between the PSD and any of the IMs 
under consideration, the Spearman rank correlation coefficient (Maritz, 1995), ps, 
for nonlinear correlations is used here. For u = 0.1, the computed values of o; are 
0.38 (0.38) and 0.80 (0.85) for the correlations between usxo and agxo (Usro and ago), 
and usso and Vero (Usro and vero), respectively; and for u = 0.2, these values are 0.56 
(0.60) and 0.69 (0.73), respectively. The efficiency of PGV as an IM is relatively 
high, and it improves with increased sliding extent, which is appreciated because 
design is mainly concerned with sliding displacements that are sufficiently large and 
may exceed the preset threshold. Nevertheless, the variability of the PSD at a given 
PGV is still considerable in comparison with the peak superstructure response of SB 
structures presented in Chap. 4. This relatively large variability is primarily attributed 
to the following reasons, which have been pointed out by Jampole et al. (2020): 


(1) The initiation of sliding is dominated by the acceleration quantities; thus, initi- 
ating sliding through a pulse with larger PGV may not be easier as the larger 
PGV could be the result of longer duration instead of a larger acceleration 
amplitude. 

(2) Although the incremental velocity of a pulse can effectively describe the sliding 
excursion resulting from a velocity pulse, the incremental velocity of the largest 
pulse may not necessarily be in close proximity to the PGV of a ground motion 
record because the value of PGV is also affected by the initial conditions 
preceding the largest pulse. 

(3) The PSD obtained from seismic excitation of a SB structure is an accumulative 
result of multiple sliding excursions initiated by large velocity pulses, especially 
for lower friction levels. The efficiency of PGV is further decreased due to the 
accumulative effect, as it is typically associated with a dominant pulse. 


Jampole et al. (2020) suggested a new IM, named EIGV, that is more effective in 
forecasting PSDs of rigid bodies exposed to earthquake ground motions. However, 
because the correlation between the PSD and PGV is acceptable and PGV is simple 
and well accepted by the engineering community as a ground motion IM, PGV is 
adopted herein. 


6.2 Critical Parameters and Their Ranges 


The parametric study presented in Chap. 5 indicates that the story number N generally 
does not influence the distribution of the equivalent lateral forces, which implies that 
the sliding displacement of a multistory SB structure should be close to that of the 
corresponding single-story SB structure with the same mass ratio and fundamental 
period. To verify this inference, the responses of a three-story and a single-story 
SB structure with the mass ratio œ = 0.75 and the fundamental periods in the x 
and y directions equal to 0.4 s were computed. Figure 6.2 compares the probability 
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Fig. 6.2 Comparison of the probability distributions of the normalized PSDs of three-story and 
single-story SB structures (o = 0.75 and T, = T, = 0.4 s): a vgxo/u = 2 m/s; b Vexo/ = 6 m/s; 


€ Vgxo/ u = 12 m/s; d vgro/u = 2 m/s; e vgro/u = 6 m/s; and f vgro/u = 12 m/s 


distributions of the PSDs normalized by the friction coefficient u of these two struc- 
tures. As can be seen, the probability distributions of the normalized PSDs of the 
three-story and single-story structures almost coincide at a given level of the PGV 
normalized by u. Therefore, single-story SB structures can be used to evaluate the 


PSDs of general SB structures. 


As already presented in Chap. 4, the governing equations of single-story SB 


structures are as follows: For the stick phases, 


T s 2 = 
lij, + 2E,@ Uy + WMUrxy = —Ugx 


s : 2 B 
li, + 2&0, u,, + WyUry = —Ugy 


The precondition for the stick phases is 


V (oi. i "Ma + (oii, + ir < (8 + tig.) Ms 
For the sliding phases, 


.. .. B 2 uu .. 

Usy + py + 25,0, + Wy.Urxy = —Ugx 

P ls is Š " 

Usx + FF (8 + Ugz) M+ QUrx = —Ugx 
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ls, 


Juss + ligc) M + Qü,y = —ü gy 


(6.4) 


(6.5) 


(6.6) 
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As can be seen from Eqs. (6.4) to (6.6), the responses of SB structures are influ- 
enced by the friction coefficient (the dynamic and static friction coefficients are 
assumed to be the same), which is a critical parameter. By dividing both sides of 
Eqs. (6.4)-(6.6) by u, it can be observed that the effects of u can be incorporated into 
both the response quantities and ground motion IM, i.e., the normalized displace- 
ment quantities, u4,(f) = uj, (f)/u, uj, (t) = Uyy(t)/M, Usx(t) = usx(t)/u, and 
Us, (t) = usy(t)/u, are dependent on the normalized IMs, vgxo/4 and vgyo/u (where 
Vgyo and vgy are the PGVs in the x and y directions, respectively), and jz is not an 
independent variable anymore. Thus, to simplify the estimation of PSDs related to 
different levels of ground motion intensity associated with various levels of u, we 
can evaluate the normalized PSDs (usx0/4, Usyo/ H, and uso/u) at various normalized 
PGVS (Vgx0/ LL, Vgyo/ H, and vVgo/ u). Equations (6.4)-(6.6) shows that cx, wy, Ex, Ey, 
a, Wg, (f), and the horizontal ground motion waveform are other parameters that may 
affect the normalized PSDs. 

The common ranges of the natural periods of the superstructure (T, = 2zt/o, 
and T, = 2zt/oy), the damping ratios (€, and &,), and the mass ratio (œ) have been 
presented in Sect. 4.2, and thus are not repeated here. The friction coefficients of 
the sliding interfaces investigated for SB structures (Barbagallo et al., 2017; Hasani, 
1996; Jampole et al., 2016; Nanda et al., 2012; Qamaruddin et al., 1986; Yegian 
et al., 2004) range from 0.07 to 0.41. Apart from very few near-fault records from 
high-magnitude earthquakes, most of the ground motions recorded have PGVs below 
1.2 m/s. On these bases, the normalized PGV is limited to 18 m/s, which is equiv- 
alent to PGV = 1.26 m/s when u = 0.07. When the normalized PGV is 1 m/s, the 
PSDs associated with jz lying in the common range are well below 0.1 m, a value 
that can be considered as the lower bound of the sliding displacement threshold in 
practice. Therefore, analyzing the cases with normalized PGV below 1 m/s is not 
necessary from a design perspective. In the following parametric study, eleven levels 
of normalized PGVs, namely 1, 1.5, 2, 4, 6, 8, 10 12, 14, 16, and 18 m/s, are consid- 
ered. These levels of normalized PGVs were achieved by adjusting the value of u 
with the ground motion records unscaled. 


6.3 Parametric Study for the Normalized Peak Sliding 
Displacements 


6.3.1 Comparison of the Responses in the Two Orthogonal 
Directions 


The coupling of the friction forces in the two orthogonal directions [i.e., their resultant 
is equal to (m + mp) (2 +ü "m causing the ground motion in one direction tends 
to decrease the friction force component in the orthogonal direction, resulting in an 
increase in the sliding displacement in that direction. However, this effect is recip- 
rocal. Assuming using enough number of ground motions, the relationship between 
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the average normalized PSD and normalized PGV in one direction that was obtained 
should be the same as that in the orthogonal direction under the circumstance of 
three-component seismic excitation. To confirm this inference, Fig. 6.3 displays the 
mean values of Us,o/4 at each level of v¿xo/u considered and the data points, (vgy0/ H, 
Usyo/{L), corresponding to the response in the y direction. These data were obtained 
from response history analyses of SB structures with T, = T, = 0.4 s and a = 0.7 
using the 180 non-pulse-like ground motion records. Figure 6.3 presents only the 
data points with 1 m/s < vao/u x 18 m/s, in order to maintain consistency with 
the range of the normalized ground motion IM considered in the x direction. It was 
found that a quadratic polynomial curve can well represent the relationship between 
mean 45,9/1L and v,,o/u; therefore, a regression curve, obtained through the use of 
a quadratic polynomial equation for fitting the data points, is displayed in Fig. 6.3. 
As can be seen in this figure, the regression curve for the relationship between the 
mean usyo/u and vgyo/u agrees well with the curve of the mean usxo/u versus Vgxo/ LL. 
Therefore, it can be inferred that the relationship between the mean normalized PSD 
and normalized PGV is essentially identical for the two orthogonal horizontal direc- 
tions. However, it is important to note that the PSDs may vary greatly between the 
two orthogonal directions for an individual ground motion, despite both directions 
having the same PGVs. In terms of statistical results, the outcomes achieved for the 
x-direction through a considerable number of ground motions can be extended to 
the y-direction. For this reason, only the response in the x direction is investigated 
hereafter. 
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Fig. 6.3 Comparison of the relationship between the normalized PSD and normalized PGV for the 
two orthogonal horizontal directions 
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6.3.2 Probability Distribution of the Normalized PSD 
at a Given Level of Normalized Ground Motion 
Intensity 


Figure 6.4 depicts the cumulative probability distribution of the normalized PSD 
(Usxo/f4 and uso/p) at four distinct levels of normalized PGV (v4,9/p. and vgo/j.). 
The data used to determine these empirical cumulative distributions were derived 
from response history analyses using the 180 non-pulse-like ground motion records 
with the structural parameters T, = T, = 0.4 s and a = 0.7. The figure clearly 
depicts that the empirical distributions are asymmetrical around the sample median 
and have lengthier tails moving towards upper values. The lognormal distribution, 
which has been extensively utilized in seismic performance assessment of building 
structures (e.g., Ruiz-Garcia & Miranda, 2007; Zareian & Krawinkler, 2007), also 
presents such type of feature and, thus, could be suitable for modeling the probability 
distributions of usxo/4 at a given level of v4,9/1 and uo/u at a given level of v,9/ 
u. The sample geometric mean and sample logarithmic standard deviation (Ang & 
Tang, 2006) are typically used to estimate the two parameters (i.e., the median and 
logarithmic standard deviation) of the fitted lognormal distribution function. For this 
study, the equations for estimating the parameters can be written as 


(usso/ i), = oo incr) (6.7a) 
zi 
(Usi0/ Mm = exp » ite) (6.76) 
zi 
— i » EO - In@sx0/Wm]. mm 
"m i Xia Bou i InWs:0/Hm]. (s 


where (usxo/U)m and (usto/U)m are the medians of usxo/u and uso/u, respectively; 
OIn(u,,o/p) Md Otn(u,,o/u) are the lognormal standard deviations of u;,9/u and usol, 
respectively; (u,,9/1.); and (uso/4); are the observed value; and n is the sample size. 
However, for v,,9/1t (and v49/u) < 2 m/s, some observed values of usxo/u (and Usio/ 
u) are 0 or very close to 0, which makes Eqs. (6.7) and (6.8) invalid because the 
natural logarithm of zero does not exist and the value computed by Eq. (6.7) will be 
dominated by the natural logarithm of a value that is near 0. Thus, for these cases, 
(Usxo/{4)m and (Us:9/[4)m are taken as the counted medians, and oj, o/u) and oj, /1) 
are estimated by Eq. (6.9) based on the assumption that the data are sampled from 
lognormal distributions. 
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Fig. 6.4 Empirical and fitted lognormal probability distributions of the normalized PSD at given 
levels of corresponding normalized PGV: a vgxo/u = 1 m/s; b vgxo/u = 2 m/s; € vgxo/u = 6 m/s; 
d vg;o/u = 12 m/s; e vgro/u = 1 m/s; f vgro/u = 2 m/s; g vgro/u = 6 m/s; h vgro/u = 12 m/s 


dang eda | EBD (6.92) 
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where (usxo/ u )s0% and (usxo/ H )gs4% are the counted median and counted 84th percentile 
Of usxo/ u, respectively; and (usro/u)so% and (usro/u)g4% are the counted median and 
counted 84th percentile of usro/ u, respectively. 

Figure 6.4 also displays the fitted lognormal distribution functions for each 
of the four levels of v,xọ/u and vgo/u. In general, the fitted lognormal distribu- 
tion agrees fairly well with the corresponding empirical distribution. The well- 
known Kolmogorov-Smirnov (K-S) goodness-of-fit tests (Ang & Tang, 2006) were 
conducted to further verify the adequacy of the lognormal distribution. Figure 6.4 
depicts the graphical representations of the K-S test with a 5% significance level. 
The figure displays that all data points for v4,9/1t (and Ve9/4) = 2, 6, and 12 m/s, are 
within the limits of acceptability (i.e., the two dotted lines in Fig. 6.4), indicating that 
the assumed lognormal distribution is acceptable. For v4;o0/14 (and vgo/u) = 1 m/s, 
due to the presence of several null values, certain points at the lower tail fall outside 
the acceptable limits; however, the practical sliding displacement threshold is much 
higher than the PSDs associated with this lower tail, thus the utility of the lognormal 
distribution remains unaffected. Concluding from the aforementioned discussions, it 
is evident that the lognormal distribution is appropriate for modeling the probability 
distributions of u;,o/u at a given level of v4,9/1,, and uso/u at a given level of vgo/ u. 
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Fig. 6.5 Comparison of the probability distributions of the normalized PSD under two- and three- 
component excitations: a Vgxo/ u = 1 m/s; b vgxo/u = 2 m/s; € vgxo/u = 6 m/s; d vgxo/u = 12 m/s; 
€ Vgro/ u = 1 m/s; f vgro/ u = 2 m/s; g Vgro/ u = 6 m/s; h vgo/u = 12 m/s 


6.3.3 Effect of the Vertical Ground Motion Component 


To investigate the influence of the vertical component on the PSD, the computation 
was carried out on the responses of SB structures exposed solely to the two horizontal 
components of the 180 non-pulse-like ground motions. Figure 6.5 shows a compar- 
ison between the normalized PSD probability distribution under the two-component 
excitation and that of corresponding three-component excitation. The normalized 
PSD probability distributions for both cases in this figure are almost the same at a 
given normalized PGV level, with only slight differences in the upper portions of 
the cumulative distribution curves. This result indicates that the effect of the vertical 
component on the PSD is negligible. Shao and Tung (1999) and Konstantinidis and 
Nikfar (2015) have also arrived at comparable conclusions regarding the sliding 
behavior of rigid bodies. 


6.3.4 Effects of the Superstructure Natural Period and Mass 
Ratio 


Figures 6.6 and 6.7 present the relationships between (Us,0/4)m and vgxo/ p, and (usio/ 
/4)m and vgo/u, respectively, for different values of Ty and a, which were determined 
by using the 180 non-pulse-like ground motion records and assuming T, = Ty. As 
shown in these figures, the trend of (u,.)/u),, changing as vgxo/u increases closely 
resembles the trend of (us;o/{4)m changing as Vg;9/j4 increases. In comparison with 
the normalized PGV, the influence of T, and œ on (u,o/u), and (Uso/{4)m is not 
so significant. To further investigate the combined effects of T, and o on (usxo/U)m 
and (usto/U)m, (usxo/ 1) and (Uso/[4)m are plotted against Ty and o in Figs. 6.8 and 
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6.9, respectively, for four representative levels of normalized PGV. This figure also 
presents the results of rigid bodies, which correspond to Ty = 0. As can be seen in 
Figs. 6.8 and 6.9, when T, < 0.4 s, the mass ratio basically has a negligible effect 
on (Usxo/[L)m and (uso/H)m; when T, > 0.4 s, the influence of a becomes slightly 
more significant, and (u,,0/1.),, and (Us0/[4)m generally increase as o increases. This 
phenomenon cannot be simply interpreted using the governing equations presented 
previously; additionally, since a larger value of o does not always lead to a larger (u,,0/ 
[L)m Of (Usto/ ),, as presented in Figs. 6.8 and 6.9, the inherent characteristics of the 
ground motion time history should have played a significant role in this general trend. 
For a given mass ratio, the values of (u;;9/1),, and (usto/U)m generally first increase 
and then decrease as T, increases, and the differences between the maximum and 
minimum values of (u,,0/14), and (Usto/H)m for Ty within the range considered range 
from 0.04 to 0.33 m and 0.02 to 0.27 m, respectively, and generally increase as the 
corresponding normalized PGV increases. From this result, we know that the PSDs 
of actual SB structures may be underestimated by relying solely on the response of 
rigid bodies. For simplicity, it is reasonable to use the maximum values of (Usx0/ 1.) 
and (Us10/4)m for the range of T, considered to conservatively estimate the PSDs of 
possible SB structures. 

Figures 6.10 and 6.11 present the relationships between Gintu, o/u) and Vgxo/H, 
and oj,(,,,/4) and Vgro/ H, respectively, for different values of T, and œ. For some 
cases when v¿xo/u (and vgro/ u) = 1 m/s, the values of (u;,9/1)50«, [and (usio/1:)50] 
are 0 or very close to 0; thus, the values obtained from using Eq. (6.9) to compute 
CIn(u,,o/ i) LANA 615,74) ] are infinite or unreasonably large. For this reason, the results 
corresponding to v4,o/u (and vgo/u) = 1 m/s are not presented in Fig. 6.10 (and 
Fig. 6.11). As shown in Fig. 6.10 (and Fig. 6.11), Cintu, o/u) [and Oinw,,o/u)] generally 
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Fig. 6.7 Relationships between (us;0//4)m and vgo/u for different values of Ty 
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Fig. 6.9 Combined effects of Ty and o on (usto/H)m: a Vgro/u = 2 m/s; b vgro/u = 6 m/s; € vgro/ u 
= 10 m/s; d vgo/u = 14 m/s 


lie between 0.4 and 0.6 except for some cases when vgxo/u (and vgo/u) = 1.5 and 
2 m/s. When the normalized PGV is small, sliding is not predominant, and the 
ground acceleration may play a more significant role than the ground velocity as the 
acceleration quantities dominate the initiation of sliding as revealed by Eq. (6.5); 
furthermore, sliding is less likely to occur for smaller values of o and larger values 
of T,, as demonstrated in Chap. 4. This explains why the values of oint, o/u) [and 
Oin(ugo/u)] are generally larger for vgxo/ 4 (and vgo/ u) = 1.5 and 2 m/s and even larger 
values are obtained when T, = 1 s anda < 0.7. 

The combined effects of T, and œ on oj(,,,/; and oj(,,,/,) are plotted in 
Figs. 6.12 and 6.13, respectively, for four representative levels of normalized PGV. 
In general, the influence of o on Ojn(u,.5/2) [and Oin(u 9/4) ] is small except for vgxo/ u 
(and vgro/ u) < 2 m/s. For any given level of v,;o/1 (and vgo/u), the maximum value 
Of e1,(,,.,/u) [and Oin(u,.9/)] is obtained at T, = 1 s; this value is 0.87 (and 1.12) 
for Vgxo/ (and vgro/ u) = 2 m/s and is around 0.63 (and 0.58) for all other levels of 
Vgxo/ u (and Vgio/ u) = 4 m/s. The value of T, at which the minimum or; [and 
€1n(u,,)/u)] is obtained varies for different levels of vgxo/u (and vgro/ u). For a given 
level of v¿xo/u (and vgo/u), the minimum value of Omu, o/u) [and om] ranges 
from 0.41 to 0.57 (and 0.38—0.53). The average value of intu, o/u) [and Oin(ugo/p)] is 
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Fig. 6.10 Relationships between Oin(u,,ọ/u) and Vgxo/ for different values of T, and a: a œ = 0.5; 
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Fig. 6.11 Relationships between ojy(u,,9/) and vgro/p for different values of Ty and o: a a = 0.5; 


bo =0.7;ca=0.9;d7,=0.1s;e Ty = 0.4 s; f T, = 0.7 s 


0.60 (and 0.58) for Vexo/ (and vgro/u) = 2 m/s and ranges from 0.49 to 0.58 (and 
0.45-0.55) for Vero/ (and vgo/ H) = 4 m/s. 

The equality 7, = T, is employed in all of the aforementioned analyses. To 
investigate the possible effect of T/T, on the PSD, the values of (usxo/4)m and (Usio/ 
/4)m corresponding to different values of T/T, are compared in Figs. 6.14 and 6.15, 
respectively, and those of omu, o/u) and Oin(u,,/.) are compared in Figs. 6.16 and 6.17, 
respectively. For the data presented in these figures, o is taken as 0.7. These figures 
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Fig. 6.12. Combined effects of T, and o on oj, 1): a Vgxo/ U = 2 m/s; b veyo/ u = 6 m/s; € Vgxo/ 
u = 10 m/s; d vgxo/u = 14 m/s 


make it clear that T,/T, has a negligible influence. Therefore, the results obtained 
for T, = T, can represent those of the other T',/T, within the range considered. 


6.3.5 Effect of Near-Fault Pulses 


Distinct pulses in near-fault ground motions affected by forward directivity may result 
in distinct sliding response characteristics as compared to ordinary ground motions. 
To investigate this effect, the (u,,9/1.),, versus Vgxo/u and (usio/p), versus vgro/ 
curves obtained using the 60 near-fault pulse-like records and the 180 non-pulse-like 
records are compared in Figs. 6.18 and 6.19, respectively. When vṣxo/u < 4 m/s (and 
Vgo/ u < 6 m/s), the values of (usxo/U)m [and (usto/U)m] corresponding to the pulse- 
like records are close to those corresponding to the non-pulse-like records. When 
Vexro/ exceeds 6 m/s (and vgo/u exceeds 8 m/s), the value of (usxo/U)m [and (us;o/ 
U)m] for the pulse-like records starts to exceed the corresponding value for the non- 
pulse-like records, and the difference increases monotonically as v,,o/u (and vgro/ 4) 
increases. To interpret the underlying reason for this phenomenon, Fig. 6.20 (T, = 
T, = 0.4 s and a = 0.7 are adopted) illustrates the ground acceleration, velocity, and 
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Fig. 6.17 Effects of T,/Ty on o1,(4,,/4) (€ = 0.7): a T; = 0.4 s; b T, —20.7s; cT, —1s 


sliding displacement time histories corresponding to the counted median of u;,o/u in 
each (non-pulse-like or pulse-like) group. By comparing the time histories presented 
in Fig. 6.20b for v4,9/u = 10 m/s, we can find that the prominent long-period velocity 
pulse in the pulse-like ground motion is the cause of the larger value of (u,,0/1),, in 
comparison with the non-pulse-like ground motion. However, when v4,o9/u is small, 
as illustrated in Fig. 6.20a for v4,9/14 = 2 m/s, the contribution of the long-period 
velocity pulse is not so significant. The simplified equation (Eq. 6.10) proposed by 
Jampole et al. (2018) for predicting the PSD of a rigid block subjected to a half- 
sine pulse which was derived from simplification of the corresponding closed-form 
solution can provide an approximate interpretation of this result. 


272 
_ 47, 1 12, ly 
Us,max = 4ug 740!» + ge Hg (6.10) 


where Us max is the PSD of the rigid block; and a, and T, are the peak acceleration 
and duration of the half-sine pulse, respectively. Dividing both sides of Eq. (6.10) 
by u and replacing a,T,, with 7 v4;/2 (where vp; is the peak velocity of the half-sine 
pulse), lead to 
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s,max i i 1 
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Fig. 6.18 Effects of near-fault pulses on (Usx0/[4)m: a œ = 0.5; b œ = 0.7; and c a = 0.9 
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Fig. 6.19 Effects of near-fault pulses on (Usi0/[)m: a æ = 0.5; ba = 0.7; and c œ = 0.9 


Equation (6.11) indicates that the quadratic relationship between normalized PSD 
and normalized PGV remains when a half-sine pulse is applied for excitation, that 
is, as the normalized PGV increases, the rate of increase of the normalized PSD with 
respect to it also increases, i.e., the effect of the prominent long-period velocity pulse 
is more significant when normalized PGV levels are higher. 

Figures 6.21 and 6.22 compare the oj,(,,.,/,) Versus Vexo/f and 0j,(,,,/,) Versus 
Vgro/ u curves, respectively, of the pulse-like records with those of the non-pulse-like 
records. As shown in this figure, when v,,o/u. (and vgo/u) x 4 m/s, the values of 
Oin(ugo/2) [aNd Otn(u,,./.)] for the pulse-like records are generally larger than those for 
the non-pulse-like records; when v4,9/1 (and vgro/ H) = 6 m/s, the value of oint, o/u) 
[and ojn(w,,o/)] for the pulse-like records does not change much and is slightly smaller 
than the corresponding value for the non-pulse-like records. Since the computed 
dispersion is partly influenced by the selected ground motion records, and typically 
there are minimal differences in the computed Gintu, o/u) [and oj,,/,)] of the two 


ground motion types, it is reasonable to expect a similar level of inherent dispersion 
for the two types of ground motions. 
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Fig. 6.21 Effects of near-fault pulses on o15(,,, 9/4): a € = 0.5; ba = 0.7; ca = 0.9 
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Fig. 6.22 Effects of near-fault pulses on o1(u,,,/4): a € = 0.5; ba = 0.7; ca = 0.9 


6.4 Fragility Curves 


From the investigations in the preceding section, we know that the influence of Tx 
and a on (u5,9/1.),, [and (usto/H)m] is limited in comparison with that of vg.o/ (and 
Vgro/ H). Therefore, in the design of SB structures, it is advisable to use the maximum 
values of (Usxo/U)m and (Uso/4)m conservatively for the common range of T,. The 
maximum (Usx0/H)m versus vg,o/t, and maximum (Us;9/{L)m Versus Vgro/ u curves are 
plotted in Fig. 6.23 for different values of a and for both the non-pulse-like and 
pulse-like ground motions. Since the curves corresponding to different values of o 
are very close to each other, for simplicity, equations for design can be developed 
solely based on the findings of æ = 0.9, which are generally larger than those of other 
values of a. It is found that a quadratic polynomial curve can well fit the relationship 
between the maximum (Us,0//4)m and Vpxo/fZ, as well as the maximum (Us;9/{4) and 
Vgro/ u, and the obtained regression formulae are as follows: 


(1) For the non-pulse-like ground motions, 


(s.o/ 1), = 0.0052(vgx0/ U)" + 0.261 (vexr0/14) — 0.2540  (6.12a) 


(Us10/ Lm = 0.0047 (vg0/ u4)" + 0.283(vg10/u) — 0.308 > 0 — (6.12b) 
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Fig. 6.23 Comparison of the design equations and numerical results for the relationships between 
the median normalized PSD and normalized PGV: a non-pulse-like ground motions; b pulse-like 
ground motions 


(2) For the near-fault pulse-like ground motions, 


(sco / 1), = 0.017 (vgx0/ U)” + 0.257(vgx0/ u) — 0.3410 — (6.13a) 


(Usto/U)m = 0.020(vei0/H)” + 0.190(vg:0/) — 0.282 > 0 (6.13b) 


where (usxo/U)m and (u,o/u), are in m, and Vpgxo/f4 and vgo/u are in m/s. 
According to the findings presented previously, replacing the subscript letter 
“x” with “y” enables the application of Eqs. (6.122) and (6.132) to the response 
in the y direction as well. As shown in Fig. 6.23, Eqs. (6.12) and (6.13) can 
well predict the corresponding relationships between the median normalized 
PSD and normalized PGV, and the coefficients of determination, R?, of these 
equations are all larger than 0.99. Further comparison of the curves determined 
by Eqs. (6.12a) and (6.12b) [and Eqs. (6.13a) and (6.13b)], as presented in 
Fig. 6.23, indicates that the relationship between the median normalized PSD 
and normalized PGV in each principal direction is close to that with respect 
to the origin. Since the median normalized PSD versus normalized PGV curve 
corresponding to Eq. (6.12b) [and Eq. (6.13a)] is slightly above that corre- 
sponding to Eq. (6.12a) [and Eq. (6.13b)], Eqs. (6.12b) and (6.132) can be used 
conservatively as unified equations for predicting the response in each principal 
direction as well as the maximum response over all the directions. 
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Itis worth mentioning here that Ryan and Chopra (2004) also proposed a design 
equation for calculating the median peak displacements, (Us:9)n, of friction pendulum 
isolators: 

TNR 4.36 T4 p(—0.99-0.201n 9) 


Ari max (vxo, Vey0) (6.14) 


where T is the isolation period, and 7 is defined as 


jc ve (6.15) 


c MAX(Vgx0, vgyo) 


wa in Eq. (6.15) is the frequency marking the transition from the velocity-sensitive to 
the displacement-sensitive region of the median spectrum of the stronger horizontal 
ground-motion components. Note that in Eq. (6.14), the PSD with respect to the 
origin is taken as the response quantity of interest, while the PGV of the stronger 
component is taken as the ground motion IM, which sets it apart from the treatment 
in this study. The comparison of the median normalized PSD versus normalized PGV 
curve, determined by Eq. (6.14) [wg = 3.05 is adopted, as done by Ryan and Chopra 
(2004), and T; is taken as 10 s such that the corresponding radius of the FP isolator 
is sufficiently large to yield the same response as that of a flat sliding system], with 
those determined by Eqs. (6.12) and (6.13) in Fig. 6.23, is presented. As can be seen, 
the curve determined by Eq. (6.14) is close to those determined by Eq. (6.13), which 
is proposed for near-fault pulse-like ground motions. This is because the 20 ground 
motions used in the response history analyses conducted by Ryan and Chopra (2004) 
were from large-magnitude earthquakes and recorded at sites near fault ruptures, the 
characteristics of which are close to the near-fault pulse-like ground motions used in 
the present study. 

As demonstrated in the preceding section, there does not exhibit a clear trend for 
the influence of the structural parameters, Tx, T ,/T ,, and o on the logarithmic stan- 
dard deviations, oj, o/u) and oj, ,/,). and the values of oj, o/u) [and oj, y] 
are generally between 0.4 and 0.6 except for some cases when vgxo/u (and vgo9/ 
4) is small. Based on these results, adopting a constant value for oj,(,,,/,) [and 
€1n(u,,/u)] 1$ reasonable in design. This value is taken as 0.55 for both oq, 7, 
and c;,(,,,/,) here, which is approximately the average of all the results computed 
using the non-pulse-like ground motion records when v,,9/u (and vgo/ u) = 2 m/s. 
As previously discussed, the dispersion for the pulse-like and non-pulse-like ground 
motions is expected to be the same; thus, the values of oj,(,,.,/,) and Oin(w,.o/n) for 
the pulse-like ground motions are also taken as 0.55. The aforementioned dispersion 
is a result of the random nature of the ground motion, which belongs to the aleatory 
uncertainty. Other sources of variability are referred to as the epistemic uncertainty, 
which is related to the lack of knowledge about the real structural properties and 
modeling approximations. Simultaneous consideration of both types of uncertainty 
involves an elaborate Monte Carlo simulation with appropriate distribution functions 
for the structural properties, which requires considerable effort. For simplicity, an 
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approximate method based on the assumption that the effects of aleatory and epis- 
temic sources are independent (FEMA, 2009) is adopted here. Assuming a recom- 
mended epistemic dispersion of 0.35 for average modeling quality, as suggested 
by FEMA P-58-1 (FEMA, 2018), the total dispersion of the normalized PSD is 
0.55? + 0.35? = 0.65. 

The fragility curve is an effective approach in assessing the seismic vulnerability 
of SB structures caused by excessive sliding, which presents the probabilities of 
exceeding a specified sliding displacement threshold at various levels of ground 
motion intensity. Since the normalized PSD (usx0/4, Usyo/ 4, and uso/ 4) at a given level 
of corresponding normalized PGV (vgx0/ 4, Veyo/ 4, and vgro/ H) follows the lognormal 
distribution, the probability, Pr, of exceeding the sliding displacement threshold, 
Ulim, for given values of PGV = pgv and jz = jio can be computed by 


Py = P(Us0 > uig| PGV = pgv, u = uo) 
= P[(U;o/Ho) > (Utim/ Mo) (PGV/p) = (pgv/uo)] 
E s (Serm = d 


(6.16) 


where U.o represents the PSD of interest; ® is the standard normal cumulative 
distribution function; the median normalized PSD, (uso/Ho)m, is computed using 
Eq. (6.12) or Eq. (6.13); and the total dispersion, Brot, is taken as 0.65, as discussed 
previously. Figure 6.24 shows the fragility curves for some typical values of u and ttim 
[Eqs. (6.12b) and (6.13a) were used in the computation], which clearly demonstrate 
the variation of Py as the PGV increases and the effects of the primary parameters. 
As shown in Fig. 6.24, for small values of uj (e.g., “im = 0.1 m) or large values of 
u (e.g., u = 0.4), there is no significant difference between the fragility curves of the 
non-pulse-like and pulse-like ground motions; as wim increases or u decreases, this 
difference becomes more significant and the SB structures subjected to pulse-like 
ground motions are more vulnerable in comparison with those subjected to non- 
pulse-like ground motions. This outcome agrees with the differences in the median 
normalized PSD and normalized PGV relationships depicted in Fig. 6.18 for both 
ground motion types. 


6.5 Conclusions 


This chapter presents a comprehensive study on the peak sliding displacements of 
SB structures subjected to three-component earthquake excitations. The PSDs in 
both the two main directions and with respect to the origin are taken into account. 
PGV is chosen as the ground motion IM because it exhibits a higher correlation with 
PSD compared to PGA and its attenuation relationship is conveniently accessible for 
design use. 
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Fig. 6.24 Fragility curves: a u = 0.1; b u = 0.2; c u = 0.3; d u = 0.4 


It is possible for an individual ground motion to exhibit significant differences in 
its PSDs between the two orthogonal horizontal directions, despite both directions 
having the same PGVs; however, on average, the relationship between the normalized 
PSD and normalized PGV is essentially identical for the two orthogonal directions. 
The effect of the vertical ground motion component on the PSD is negligible. The 
probability distributions of us,9/ at a given level of v4,9/p and uso/u at a given level 
Of vgo/u can be modeled by the lognormal distribution. The relationship between 
(Usxo/{4)m and vg,o9/p. and that between (usro/U)m and vg;o/u are close to each other. 
The influence of T,, T,/T,, and o on (u;;o/1), and (Uso/4)m is insignificant; thus, 
it is appropriate to conservatively use the maximum values of (u,,5/14),, and (ustol 
/4)m for the common ranges of T,, T,, and o in the design of SB structures. The 
lognormal standard deviation, Onu, o/u) [and Oinu,.o/u)], generally lies between 0.4 
and 0.6 except for some cases when v4,o0/1. (and vgo/ u) is below 2. 

When the normalized PGV is small, the values of (us.0//4)m [and (u;;0/1)5] corre- 
sponding to the pulse-like records are close to those corresponding to the non-pulse- 
like records. When the normalized PGV exceeds a certain value (approximately 


100 6 Peak Sliding Displacements of Sliding Base Structures Under ... 


6-8 m/s), the value of (Usx0/{4)m [and (Usi0/4)m] for the pulse-like records starts to 
exceed the corresponding value for the non-pulse-like records, and the difference 
increases monotonically as v4,o/u (and vgo/ 4) increases. The difference in the value 
Of Onu, o/u) [aNd Oin(u,.o/u)] for the two types of ground motions is small. 
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